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ABSTRACT 


This  interim  report  describes  the  work  to  date  on 
development  of  the  ATR  computer  code  -  a  user-oriented 
code  for  calculating  quickly  and  simply  radiation  environ¬ 
ment  problems  at  all  altitudes  in  the  atmosphere.  The 
ATR  code  is  based  on  parametric  models  of  a  comprehen¬ 
sive  data  base  of  air  transport  results  which  was  generated 
using  rigorous  transport  techniques  for  infinite  homogeneous 
air.  The  effects  of  an  air-ground  interface  and  non-uniform 
air  density  are  treated  as  perturbation  corrections  on 
homogeneous-air  results.  The  data  base  has  been  generated 
for  neutrons,  secondary  gamma-rays,  high-energy  photons 

and  low  energy  photons  as  a  function  of  space,  energy  and 

o 

angle  out  to  ranges  of  550  gm/cm  of  air.  Parametric 
models  have  been  developed  for  all  of  this  data  base  except 
photon  angular  distribution  and  a  version  of  ATR  has  been 
completed  for  neutron  and  secondary  gamma- ray  transport. 
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I. 


INTRODUCTION 


This  report  describes  research  performed  at  Science  Applications, 
Inc.  over  the  past  year  to  develop  parametric  models  for  efficiently  cal¬ 
culating  radiation  transport  in  the  atmosphere.  The  fundamental  problem 
being  solved  is  to  represent  the  free-field  radiation  environment  for 
monoenergetic  sources  by  parametric  equations  so  that  free-fields  for 
arbitrary  source  spectra  can  be  calculated  quickly  by  simple  folding 
operations. 

Our  general  approach  was  to  develop  detailed  parametric  models 
by  curve  fits  to  rigorous  transport  data  for  infinite  homogeneous  air,  and 
to  include  the  effects  of  an  air-ground  interface  and  exponential  air  density 
as  perturbation  corrections.  The  end  result  of  this  effort  is  the  ATR 
computer  code  which  has  been  designed  as  a  user-oriented  code  to  calcu¬ 
late  a  wide  variety  of  radiation  environment  problems  at  all  altitudes  in 
the  atmosphere. 

The  work  to  date  includes:  (1)  the  development  of  a  comprehensive 
data  base  of  homogeneous  air  transport  results  as  a  function  of  source 
energy,  space,  energy  and  angle  for  neutrons,  secondary  gamma-rays, 
high-energy  photons  and  low-energy  photons;  (2)  parametric  models  for 
the  entire  data  base  except  for  photon  angular  distribution;  (3)  corrections 
to  homogeneous  air  results  for  neutrons  and  secondary  gamma-rays  and; 
(4)  a  completed  version  of  ATR  for  neutrons  and  secondary  gamma-rays. 
Detailed  discussions  of  these  items  are  given  in  subsequent  sections. 
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II.  DATA  BASE  FOR  RADIATION  TRANSPORT  IN  HOMOGENEOUS 

AIR 

Although  many  volumes  of  data  for  radiation  transport  in  the  air 
have  been  compiled  over  the  years  (see  Appendix  A)  much  of  it  has  been 
out-dated  because  of  newer  cross-section  data.  Even  more  recent  results 
are  often  available  only  in  tabulated  form  or  only  for  a  selected  range  of 
values  in  the  independent  variables.  All  of  these  factors  make  it  difficult 
to  assemble  a  consistent  comprehensive  data  base  from  existing  results 
which  is  suitable  for  production-processing  by  computer.  For  this 
reason,  we  have  generated  a  complete  and  uniform  data  base  at  SAI  using 
rigorous  transport  codes  based  on  the  most  recent  cross-section  data 
available.  These  calculations  will  be  described  in  this  section  including 
selected  comparisons  with  other  calculated  results  and  experiments.  It 
is  anticipated  that  the  complete  data  base  will  be  available  in  published 
form  at  the  end  of  this  contract  year. 

2. 1  Neutron  and  Secondary  Gamma-Rays 

The  data  base  for  neutron  and  secondary  gamma-ray  transport 
was  based  on  Straker’s  results  (^published  in  ORNL  4464.  These  data  were 
considered  to  be  the  most  definitive  and  comprehensive  set  of  calculations 
available  at  the  time  and  have  been  thoroughly  documented  including 
detailed  comparisons  with  integral  experiments.  The  SAI  data  base  was 
generated  using  the  identical  P5S16  air  cross-sections  of  ORNL  4464 
provided  by  Straker  in  22  neutron  groups  and  18  gamma  groups. 


3. 


The  SAI  calculations  were  made  with  ANISN  for  18  monoenergetic 
source  bands  and  for  a  maximum  range  of  5000  meters  at  a  uniform  air 

O 

density  of  1.11  mg/cm  .  Since  these  calculations  were  made,  recently 
evaluated  cross-sections  for  nitrogen  and  oxygen  have  become  available 
from  the  work  of  Young.  Selected  calculations  were  performed  using 
Young’s  cross-sections  to  determine  their  effect  on  the  existing  data 
base.  The  dose  versus  distance  results  for  neutron  and  secondary  gam¬ 
ma-rays  are  shown  in  Fig.  1  and  Fig.  2,  respectively,  for  a  14  MeV 
neutron  source.  The  largest  discrepancy  occurs  for  secondary  gamma- 
rays  with  the  results  obtained  using  Young's  cross-sections  being  lower 
than  Straker's  by  as  much  as  a  factor  of  2.  The  source  of  this  difference 
is  indicated  in  Fig.  3  which  shows  Young's  gamma-production  cross- 
sections  being  lower  than  Strakers  by  a  factor  of  two.  The  question  of 
which  is  the  more  accurate  cross-section  data  seems  unclear  at  this 

time  because  of  conflicting  conclusions  between  theory-experiment  com- 

(17) 

parison.  Secondary  gamma-ray  results  from  the  HENRE  experiment  are 
in  good  agreement  with  Straker's  cross-sections,  but  the  recent  re- 

/t 

analysis  of  LNg  experiments  by  Reynolds  is  in  better  agreement  with 
Young's  cross-sections.  It  is  aniicipated  that  this  discrepancy  will  be 
resolved  soon  and  allow  the  data  base  for  ATR  to  be  updated  accordingly. 

2.  2  High-Energy  Photons 

The  photon  transport  data  base  is  needed  for  a  wide  range  of 
energies  extending  from  about  10  keV  up  to  10  MeV.  Two  separate 
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Fig.  1.  4 »rRz  neutron  dose  versus  range  for  a  point  isotropic 
14  MeV  source  in  infinite  homogeneous  air. 


Fig.  2.  4*R  secondary  gamma-ray  dose  versus  range  for  a  point 
isotropic  14  MeV  source  in  infinite  homogeneous  air. 


/ 

l 

transport  calculations  were  made  in  order  to  achieve  an  optimum  group 
structure  for  both  the  high  and  low  photon  energy  ranges.  The  high- 
energy  photon  transport  data  base  was  calculated  using  DTFXRAY  in  the 
P5S16  approximation  with  an  analytical  first-collision  source.  The 
group  structure  was  the  same  18  group  representation  from  .02  to  10 
MeV  used  for  the  secondary  gamma-rays  and  the  maximum  range  was 
5000  meters  in  air  of  density  1.11  mg/cm  .  GAMLEGX  with  pair  pro¬ 
duction  was  used  to  prepare  multigroup  cross-sections.  The  air  com¬ 
position  by  weight  was  0.  2318  for  oxygen,  0.7553  for  nitrogen  and  0.0129 

(78) 

for  argon.  Input  absorption  cross-sections  were  taken  from  Hubbell 

for  energies  above  1  MeV  and  from  McMaste/79^  for  energies  below  1  MeV. 

Scattering  cross-sections  are  calculated  internally  in  GAMLEGX  by  the 
Klein-Nishina  formula.  All  photon  transport  calculations  were  made 
without  coherent  scattering  since  GAMLEGX  can  only  treat  coherent 
scattering  as  isotropic.  This  approximation  tends  to  be  poor  at  low 
energies  where  coherent  scattering  is  the  most  important. 

Before  generating  the  high  energy  photon  data  base,  several 
test  calculations  were  performed  to  test  the  overall  sensitivity  of  the 
results  to  choice  of  cross-section  data  and  transport  code.  Figure  4 
shows  a  comparison  of  dose  versus  distance  results  calculated  with 
different  codes  and  cross-sections  sets.  A  maximum  difference  of 
about  20%  occurs  at  100  gm/cm^,  with  the  ANISN,  DTFXRAY  and  OGRE 
results  being  in  very  close  agreement  with  each  other. 


8. 


•  ■  ,.vil 


X 


DISTANCE  (GM/CM2) 


Fig.  4.  Comparison  of  dose  vs  distance  calculations 
for  a  prompt  fission  source. 
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Figures  5  and  6  are  typical  computer-generated  plots  taken  from 
the  high-energy  photon  data  base  and  illustrate  the  general  transport 
behavior.  Figure  5  is  a  plot  of  the  scattered  fluence  versus  distance 
for  each  of  the  18  sources.  The  systematic  decrease  in  scattered 
fluence  near  the  maximum  range  (  ~550  gm/cm  )  is  due  to  leakage 
out  of  the  system  caused  by  the  vacuum  boundary  condition.  All  para- 
metric  fitting  was  extended  out  to  only  ~530  gm/cm  to  avoid  errors 
due  to  this  flux  depression.  Figure  6  is  a  plot  of  the  scattered  energy 
spectra  for  an  8.0-6.  5  MeV  source  at  distances  of  1,  3  5,  7,  and  10 
MFP  of  the  source  energy.  These  data  all  show  a  smal.  peak  at  0.  5 
MeV  due  to  the  annihilation  quanta  in  pair-production. 

2. 3  Low-Energy  Photons 

The  low-energy  photon  transport  data  base  covers  the  energy 
range  from  300  to  10  keV  and  was  generated  using  DTFXRAY  as 
described  in  Section  2.2.  The  McMaster  cross-section  compilation 
was  used  for  energies  below  1  MeV,  and  was  selected  as  being  the 
most  accurate  for  air  transport  based  on  a  series  of  sensitivity  calcu¬ 
lations  and  comparisons.  The  largest  differences  in  basic  cross- 
section  data  for  the  three  elemental  constituents  of  air  being  considered 
occur  in  the  absorption  values  for  nitrogen  at  energies  below  100  keV. 
Six  different  compilations  of  cross-section  data  were  compared  for  the 
energy  range  20-100  keV  and  are  given  in  Table  I.  The  maximum 
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TABLE  I. 


Photoelectric  cross-sections  in  N  from 
different  compilations.  Values  in 
barns/atoms. 


E 

(keV) 

1968 

LRL 

Plechaty 

1969 

LRL 

Me  Master 

1969 

NBS 

Hubbell 

1967 

LASL 

Storm 

1967 

LRL 

McMaster 

1957 

NBS 

Gradstein 

20 

7.30 

7.90 

8.26 

8.02 

8.70 

8.21 

30 

1.86 

2.06 

2.23 

2.06 

2.31 

2.15 

40 

0.706 

0.792 

0.878 

0.790 

0.898 

0.81 

50 

0.334 

0.379 

0.425 

0.377 

0.432 

0.38 

60 

0.182 

0.208 

0.238 

0.204 

0.238 

0.21 

80 

0.  0702 

0.0810 

0.094 

0. 0790 

0.0932 

0.082 

100 

0.0339 

0.0393 

0.0454 

0.0373 

0.0453 

_ 

0.041 

values  occur  in  1968  Plechaty  compilation  and  the  minimum  values  for  the 

1969  Hubbell  compilation.  The  relative  difference  between  the  maximum 

and  minimum  absorption  values  at  60  keV  is  approximately  30%,  which 

results  in  a  100%  difference  in  the  total  calculated  fluence  at  20  MFP.  In 

order  to  select  the  most  appropriate  cross-section  set,  comparisons  were 

(23) 

made  between  deep-penetration  photon  transport  experiments  in  LNg  and 

calculations  based  on  the  discrete  ordinates  code  DTFXRAY  and  the  Monte 

Carlo  code  OGRE.  The  experiment  consisted  of  measuring  the  forward 

241 

number  current  as  a  function  of  distance  from  an  Am  isotopic  source 
(59. 57  keV  gamma  energy)  in  a  LNg  dewar  that  was  16  feet  high  by  8  feet 
in  diameter.  The  comparison  of  results  are  shown  in  Fig.  7.  The  three 


4wR2J*(R)  (PHOTONS/SOURCE  PHOTON) 


0  3  6  9  12  15  IB  21 


OISTANCE  FROM  SOURCE  IN  MFP 

Fig.  7.  Comparison  of  different  calculations  with  LN? 
photon  transport  experiment.  £ 
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DTFXRAY  calculations  used  the  maximum  (1968-Plechaty),  minimum 
(1969-Hubbell)  and  average-value  (1969-Mc Master)  cross-section  sets. 

The  OGRE  calculations^^  also  used  the  average-value  1969-McMaster 
compilations.  Figure  7  demonstrates  good  agreement  between  experi¬ 
ment  and  theory  for  the  1969-McMaster  absorption  cross-sections.  In 
addition,  there  is  also  good  agreement  between  OGRE  and  DTFXRAY  for 

the  same  cross-sections.  The  OGRE  calculation  used  importance 
5 

sampling  with  10  histories,  and  resulted  in  statistical  errors  on  the 
forward  current  of  several  percent  at  deep  penetrations. 

Additional  comparisons  between  calculations  were  made  for  photon 
transport  in  air  and  are  illustrated  in  Fig.  8  for  the  energy  fluence  build¬ 
up  as  a  function  of  source  energy  at  10  MFP.  The  OGRE^^and  DTFXRAY 
results  agree  within  10%  for  most  energies,  but  differ  from  the  Kaman 

(19) 

Nuclear  Monte  Carlo  results  by  as  much  as  a  factor  of  2.  The  source  of 
this  large  discrepancy  is  not  known,  however,  the  small  number  of  initial 
histories  (2000)  and  lack  of  importance  sampling  in  the  Kaman  Nuclear 
calculations  suggests  that  large  statistical  errors  could  be  present  at 
deep  penetrations. 

Computer  plots  representative  of  the  low-energy  photon  data  base 
generated  with  DTFXRAY  are  shown  in  Figs.  9  and  10.  Figure  9 
illustrates  the  scattered  fluence  vs  distance  for  each  of  the  18  sources 
from  10-300  keV;  Fig.  10  shows  plots  of  scattered  photon  spectra  vs 
energy  for  a  300-280  keV  source  and  for  eight  distances  from  1  to  20  MFP. 
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ENERGY-FLUENCE  BUILD  UP 


T 


T 


T 


T 


—  —  -  SAl  (DTFXRAY) 

—  SAl  (OGRE) 

—  —  KAMAN  NUCLEAR  (MONTE  CARLO) 

—  UNITED  NUCLEAR  (MONTE  CARLO) 


20  40  60  80  100  200  300 

SOURCE  ENERGY (KEV> 


Fig.  8.  Comparison  of  build  up  factor  calculations  vs  source 
energy  for  low  energy  photons. 
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III.  PARAMETRIC  MODELS  OF  THE  DATA  BASE 

The  fundamental  problem  being  solved  is  to  represent  the  free- 

field  radiation  environment  <£(E  .  E,Q,  r)  for  monoenergetic  sources 

s 

E  by  parametric  equations  so  that  free-fields  for  an  arbitrary  source 
s 

spectrum  W(ES)  can  be  generated  simply  and  quickly  by  the  folding 
integral: 

$(E,  r,  O)  =  /w(E)  <p(E.  E,  R,  r)  dEo 
y  »  s  s 

Folding  monoenergetic  source  data  from  raw  flux  dump  tapes 

is  used  frequently  to  calculate  radiation  transport,  however,  computer 

storage  is  quickly  exausted  in  most  systems  before  a  complete  data  base 

can  be  stored.  The  neutron,  secondary  gamma-ray  and  photon  data  base 

7 

being  used  at  SAI  contains  nearly  10  flux  values.  In  addition  to  storage 
problems,  the  manipulation  of  tapes  is  slow  and  often  cumbersome.  The 
aim  of  the  parametric  modeling  approach  is  to  reduce  the  size  of  the  cal¬ 
culations!  problem  considerably.  This  reduction  would  allow  air  transport 
problems  to  be  calculated  from  a  self-contained,  fast  running,  easy  to  use 
computer  code.  More  emphasis  can  be  placed  on  flexibility  in  the  different 
types  of  problems  that  can  be  calculated,  in  addition  to  designing  a  code 
oriented  to  the  nerds  of  the  user. 

The  one  dominant  factor  which  makes  parametric  modeling  of 
a  large  data  base  feasible  is  the  slowly  -  changing,  asymptotic  behavior 
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characteristic  of  radiation  transport  in  large  homogeneous  systems. 
Thus,  although  structure  may  exist  in  energy  spectra,  at  a  given  range, 
the  propagation  of  this  structure  at  other  ranges  occurs  in  a  systematic 
slowly -changing  manner.  This  behavior  was  illustrated  in  the  previous 
plots  of  Section  II. 

The  general  approach  used  is  to  describe  the  free-field  quantities 
in  terms  of  shape  functions  and  the  areas  or  fluences  under  these  shapes. 
This  approach  makes  the  different  parts  of  the  problem  separable  and 
allows  greater  flexibility  in  manipulating  and  displaying  the  data  to 
determine  optimum  parametric  equations.  The  general  form  of  the 
equations  used  in  ATR  to  represent  the  angular  flux  is: 


(E  ,  r,  E, 
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R(Es,E,r,n)  ^  - - 

Jo(V^r,  E,  0)  dQ 
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if 
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t* 
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and 

o 

<p  -  angular  flux;  particles/cm  •  MeV«  steraiian*  source  neutron. 

F  -  fluencc  functions;  particles/cm  ‘source  neutron. 

S  -  Scalar  flux  shapes;  fraction/McV,  integral  over  energy  equal  unity. 

R  =  ratio  of  angular  to  scalar  flux  at  each  energy  E  for  a 

given  angular  direction,  range  and  monoc»nergetic  source 
energy;  fraction/steradiar  integral  over  solid  angle  equal  unity. 

//(E)  =  total  neutron  cross  section  (cm"*)  at  energy  E  and  air 

2 

density  of  1.11  mg/cm 

=  solid  angle  corresponding  to  the  zero-degree  direction. 


The  parametric  modeling  consisted  of  making  curve  fits  to  three  basic 
quantities  for  each  monoenergetic  source  band:  (1)  scattered  fluences 
FSCATT^’  8Cattered  spectrum  shapes  S(E,R)  and;  (3)  angular  ratio 
functions  R(E,r,  Q). 

3. 1  Fluence 

Each  set  of  fluence  data  was  fit  with  parametric  equations  for 
the  entire  range  of  180  spatial  Intervals  from  0*550  gm/cm  of  air  using 
6  coefficients  for  each  of  the  18  source  energies.  No  simple  functional 
form  was  found  to  fit  all  of  the  data  to  the  desired  accuracy  (+  2%)  for  all 
source  energies.  Two  basic  functions  were  used  for  the  fits  (in  semi- log 
space);  one  was  a  coupled  quadratic  and  the  other  the  following: 


*(R) 
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Figures  11  through  14  are  plots  of  the  raw  data  and  the  parametric  fit 
for  fluences  of  neutrons,  secondary  gamma-rays,  high  energy  photons 
and  low  energy  photons,  respectively.  This  comparison  is  typical  of 
all  fits  to  fluence  data  and  shows  typical  errors  in  the  fit  of  several 
percent. 

3.2  Energy  Distribution 

The  spectral  shape  functions  were  more  complicated  to  deter¬ 
mine  since  simultaneous  variations  in  energy  and  distance  occur  for 
each  source  energy.  Various  plots  of  the  data  were  first  studied  to 
determine  where  systematic  behavior  occured  and  what  functional  forms 
might  be  appropriate.  Figure  15,  which  illustrates  the  general  behavior, 
is  a  plot  of  normalized  neutron  spectral  data  for  ten  distances  from 

0-550  gm/cm  for  a  14  MeV  neutron  source.  The  nearly  constant 

o 

asymptotic  shape  is  apparent  for  ranges  ^100  gm/cm  .  Two- 
dimensional  fits  in  energy  and  distance  were  made  to  these  data  by  first 
fitting  the  data  with  piece-wise  continuous  functions  versus  energy  for 
each  radius,  and  then  fitting  each  resulting  coefficient  as  a  function  of 
radius.  Separate  piece-wise  continuous  functions  were  used  for  the  fits 
versus  energy  to  minimize  the  total  number  of  parameters  required. 
Thus,  only  6  coefficients  were  required  in  most  cases  to  fit  the  varia¬ 
tion  in  energy  for  a  given  distance  and  source  energy.  In  log-log  space, 
the  functions  were  (1)  a  linear  term  for  energies  below  E  £  0. 1  MeV, 
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Fig.  12.  Secondary  gamma-ray  fluence  vs  range  for 
12.2-15  MeV  neutron  source. 
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Fig.  13.  Parametric  fit  to  scattered  fludnce  for  a 
6. 5-8.0  MeV  photon  source. 
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Fig.  14,  Parametric  fit  to  the  scattered  fluence  for  a 
220-260  keV  photon  source. 
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(2)  a  quadratic  for  0.1  MeV  <  E  <  E  and  (3)  a  constant  for  E  .  The 
resulting  six  coefficients  for  each  radius  were  then  fit  vs  radius  using 
the  same  6  parameter  functions  used  in  the  fluence  fits.  Thus, 
approximately  36  total  coefficients  were  typically  used  to  describe  the 
spectral  shape  vs  distance  for  each  source  energy  giving  a  total  of 
36  x  18  648  coefficients.  The  actual  number  of  coefficients  from  these 

fits  was  usually  much  less  due  to  the  nearly  constant  spectral  shape  for 
the  last  80  or  so  radii.  The  parametric  fits  then,  reduced  the  original 
<f»(18  x  22  x  180)  array  of  71,280  neutron  fluxes  by  more  than  a  factor 
of  100  to  approximately  700  coefficients  and  several  equations. 

The  uncertainty  in  the  parametric  fits  was  small  (~10%)  except 
at  energies  where  "peaked”  structure  existed.  Errors  at  these  points 
were  usually  systematic  as  a  function  of  distance  and  were  corrected  by 
a  separate  correction  subroutine  which  fit  the  difference  between  the 
raw  data  and  parametric  representation  at  these  points.  The  features 
of  this  subroutine  allow  fine-tuning  of  the  fit  and  can  also  be  used  to 
convert  ATR  from  one  data  base  to  another  without  having  to  go  through 
the  entire  parametric  fitting  procedure.  This  in  fact  was  done  early  in 
the  program  to  convert  fits  from  a  LASL  flux  tape  to  the  existing  data 
base. 

Figures  16  through  21  illustrate  neutron  spectra  calculated  by 

the  parametric  equations  in  ATR  for  a  14  MeV  source  at  distances  from 

o 

300-4800  meters  in  air  at  a  density  of  1.11  mg/cm  .  The  solid  histogram 
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Fig.  17.  Comparison  of  neutron  spectra  calculated  by  ATR  and 
the  ANISN  data  base  for  a  12. 2-15  MeV  neutron  source 
at  600  meters  in  air  of  density  1.11  mg/cm^. 
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Fig.  18.  Comparison  of  neutron  spectra  calculated  by  ATR  and  the 
ANISN  data  base  for  a  12. 2-15  MeV  neutron  source  at 
1200  meters  in  air  of  density  1.11  mg/cm^. 
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Fiji.  19.  Comparison  of  neutron  spectra  calculated  by  ATR  and  the 
ANISN  data  base  for  a  12. 2-15  MeV  neutron  source  at 
1800  meters  in  air  of  density  1. 11  mg/cm^. 
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Fig.  20.  Comparison  of  neutron  spectra  calculated  by  ATR  and  the 
ANISN  data  base  for  a  12.2-15  MeV  neutron  source  at 
3600  meters  in  air  of  density  1. 11  rag/cm^. 
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Fig.  21.  Comparison  of  neutron  spectra  calculated  by  ATR  and 
the  ANISN  data  base  for  a  12.2-15  MeV  neutron  source 
at  4800  meters  in  air  of  density  1.11  mg/cm3. 
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is  the  raw  ANISN  data,  the  dashed  histograms  are  ATR  values  for  the 
ANISN  group  structure  and  the  smooth  curve  is  the  continuous  ATR  repre¬ 
sentation.  The  differences  between  ATR  and  ANISN  are  typically  less 
than  10%. 

Spectral  shapes  for  secondary  gamma-ray  data  displayed  a  similar 
asymptotic  behavior  vs  distance  as  neutron  spectra  and  are  illustrated  for 
a  14  MeV  neutron  source  in  Fig.  22  for  normalized  data.  Two-dimensional 
fits  to  these  shapes  were  made  using  a  fourth  order  polynomial  in  energy 
whose  coefficients  were  then  fit  vs  distance  using  one  of  the  six-parameter 
fluence  type  functions  described  in  the  previous  section.  The  structure  in 
the  shapes  at  short  ranges  was  taken  into  account  by  using  the  correction 
routine  to  fit  differences  between  the  polynomial  fit  and  raw  data.  Final 
errors  in  the  parametric  fits  of  spectral  shapes  were  less  than  10%. 

Figures  23  through  28  are  plots  of  secondary  gamma-ray  spectra  from 
ATR  and  ANISN  for  a  14  MeV  neutron  source  and  ranges  of  300  m  to  4800  m 
in  air  at  density  1.11  mg/cm  . 

Normalized  spectra  shapes  for  8  radii  are  shown  in  Figs.  29  and 
30  for  representative  high  and  low-energy  photon  sources,  respectively. 
Unlike  neutron  and  secondary  gamma -ray  transport  data,  the  photon  spectra 
displayed  a  reasonably  smooth  asymptotic  behavior  which  lends  itself  to 
a  different,  more  economical  fitting  approach.  The  procedure  used  for 
each  source,  to  model  both  high  and  low-energy  photon  transport,  was  to 
generate  spectral  shapes  for  all  distances  from  a  reference  shape  obtained 
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Fig.  22.  Secondary  gamma-ray  spectra  calculated  ey  ANISN 

for  a  12.2-15  MeV  neutron  source  at  different  ranges 
from  0-550  gm/cm^  of  air.  Data  normalized  at  low 
energies. 
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Fig.  23.  Comparison  of  secondary  gamma-ray  spectra  calculated 
by  ATR  and  the  ANISN  data  base  for  a  12. 2-15  MeV  neu¬ 
tron  source  at  300  meters  in  air  of  density  1. 11  mg/cm^. 
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Fig.  24.  Comparison  of  secondary  gamma-ray  spectra  calculated 
by  ATR  and  the  ANISN  data  base  for  a  12. 2-15  MeV  neu¬ 
tron  source  at  600  meters  in  air  of  density  1.11  mg/cm^. 
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Fig.  26.  Comparison  of  secondary  gamma- ray  spectra  calculated 
by  ATR  and  the  ANISN  data  base  for  a  12.2-13  MeV  neu¬ 
tron  source  at  1800  meters  in  air  of  density  1. 11  mg/cm3. 
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Fig.  27.  Comparison  of  secondary  gamma-ray  spectra  calculated 
by  ATR  and  the  ANISN  data  base  for  a  12.2-15  MeV  neu¬ 
tron  source  at  3600  meters  in  air  of  density  1. 11  mg/cm3 


Fig.  28.  Comparison  of  secondary  gamma-ray  spectra  calculated 
by  ATR  and  the  ANISN  data  base  for  a  12.2-15  MeV  neu¬ 
tron  source  at  4800  meters  in  air  of  density  1.1 1  mg/cm 
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Fig.  30.  Normalized  spectral  shapes  for  a  photon  source 
of  260-220  keV  at  distances  of  1  to  20  MFP. 


by  curve-fitting  the  spectrum  vs  energy  at  one  distance.  The  fixed 
distance  was  usually  chosen  to  be  at  least  several  mean  free  paths  from 
the  source. 

For  high  energy  photon  sources,  the  reference  shape  was  obtained 
by  curve-fitting  a  cubic  equation  in  log-log  space  to  the  spectral  data  for 

ENORM“E  <EMAX  '  The  fluX  at  t*lis  minimum  energy  point 
was  assumed  to  be  constant,  and  at  the  maximum  energy  point  E^^  > 
was  fit  vs  distance  using  one  of  the  two  polynomials  used  in  the  fluence 
fits.  Spectral  shapes  at  an  arbitrary  distance  were  generated  by  linear 
interpolation  from  the  reference  shape.  Thus,  in  the  region  above  0.5  MeV 
where  spreading  occurs  in  the  spectral  shapes  vs  distance  and  energy,  an 
arbitrary  spectral  shape  is  calculated  from: 


S(E,  r)  =  exp  |  In  SQ(E,  rQ)  + 
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where  Ej  =0.5  MeV  energy  point. 

Sq(E,  rQ)  =  exp  |  Aj+AglnE+AgilnE)2  +  A4(lnE)3  j 
is  the  reference  shape, 

(  2^-1 
S(emax»  =  exP  I  Bi+B2lnr+Bg(lnr)  +B4(lnr)  +Bg(lnr) 

is  the  fit  vs  distance  for  the  flux  at  E^j^  . 
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The  differences  between  the  parametric  shapes  and  the  raw  data  were 
small  except  at  l ho  0.5  McV  peak.  This  difference  was  reduced  to  less 
than  10%  for  all  ranges  except  very  close-in  by  introducing  a  constant 
correction  at  0.  5  MeV. 

A  similar  scheme  to  the  above  was  used  for  low  energy  photons 
except  for  energies  below  the  normalization  point  where  the  data  were 
not  constant  (Fig.  30)  but  exhibited  a  nearly  constant  displacement  vs 
energy  for  different  radii.  A  constant  displacement  vs  energy  was 
assumed  for  this  region  of  the  spectra  and  was  obtained  in  parametric 
form  as  a  function  of  distance  by  fitting  a  quadratic  at  one  energy  at  the 
fourth  energy  point  (  27  keV).  A  cubic  equation  also  worked  well  for 

fitting  the  reference  shape  in  the  low-energy  photon  data  base. 

An  arbitrary  spectral  shape  for  a  given  source  energy  was 
calculated  in  an  analogous  manner  to  the  above  equation  for  E  ^  ENqRM. 
For  E  <  EnqRM>  the  reference  shape  was  shifted  by  an  amount  deter¬ 
mined  from  the  fit  vs  distance  for  the  fourth  energy  point,  i.e.  for 

E  <  ENORM  ’ 


S(E,  r) 


exp  |  In 


VE’  ro> 


t  In 


S(E4,r) 


In 


S(E 


[•ro'l| 
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where  Sq(E,  rQ)  is  the  reference  shape  at  rQ  determined  by  a  cubic 
fit  vs  energy,  and  the  fourth-point  fit  vs  distance  is 


S(E4,  r)  =  exp 


Cj  +  Cg  inr  +  Cg 


(lnr)2 


3.3  Angular  Distribtuions 

The  energy-angle  coupling  of  the  scattered  flux  was  taken  into 
account  by  parametrizing  the  following  ratio  function: 

4>(ES,  r,  E,  U) 

R(Eg,  E,  r,  li)  y«<es,  r,  E, 

where  R  is  the  ratio  of  the  angular  flux  to  the  scalar  flux  at  each  energy 
E  for  a  given  angular  direction,  range  and  monoenergetic  source  energy. 
The  use  of  the  ratio  function  to  represent  angular  dependence  tends  to 
"smooth  out"  the  energy  structure  and  also  maintains  the  slowly-varying 
behavior  of  the  flux  as  a  function  of  distance  and  angle.  The  "smoothing" 
occurs  because  the  qualitative  shape  of  a  scalar  flux  spectrum  does  not 
differ  grossly  from  the  spectrum  shape  for  each  angle.  These  factors 
allowed  the  energy-angle  coupling  to  be  modeled  with  a  minimum  number 
of  parameters. 
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Parametric  modeling  for  energy-angle  coupling  has  been  com¬ 
pleted  for  m  atrons  and  secondary  gamma-rays,  and  is  discussed  below. 
The  general  behavior  of  the  energy -angular  flux  is  shown  in  Figs.  31 
through  35  at  100  meter  and  2400  meter  for  a  12.2-15  MeV  neutron 
source.  Figure  36  illustrates  the  variation  of  angular  flux  vs  distance 
for  the  6.36-8. 19  MeV  detector  group  at  6  angles.  These  plots  were 
generated  by  computer  from  the  raw  data  base  and  demonstrate  the 
variation  in  energy  spectra  for  each  of  the  17  angles  used  in  the  ANISN 
calculations.  A  significant  feature  of  these  plots  is  the  irregular 
behavior  for  small  flux  components  which  tends  to  occur  in  the  back¬ 
ward  directions  becoming  negative  for  some  cases.  This  behavior  is 
apparently  due  to  the  truncation  of  the  Legendre  expansion  at  P5  in  the 
angular  scattering  representation  of  the  discrete  ordinates  calculations. 
These  irregular  components  make  it  difficult  to  initiate  a  curve  fitting 
scheme  on  a  production  basis,  however,  they  should  be  included  for 
completeness  in  any  model.  Therefore,  as  a  first  step  toward  modeling 
the  energy-angle  coupling,  a  new  angular  flux  tape  was  created  in  which 
the  negative  fluxes  and  irregular  behavior  were  corrected  by  extrapolation 
and  interpolation.  This  procedure  was  carried  out  subject  to  the  con¬ 
straint  that  the  integrated  fluence  was  conserved  for  each  angle.  The 
important  differential  flux  data  in  the  forward  four  angles  (0-40°)  was 
essentially  unchanged  by  this  operation. 
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Fig.  31.  Neutron  angular  flux  spectra  calculated  with  ANISN 
at  17  angles  for  a  12.2-15  MeV  neutron  source  at  a 
range  of  100  meters  in  air  of  density  of  1. 11  mg/cm3. 
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Fig.  32.  Neutron  angular  flux  spectra  calculated  with  ANISN 
at  17  angles  for  a  12.2-15  MeV  neutron  source  at  a 
range  of  2400  meters  of  a»r. 
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Fig.  34.  Secondary  gamma-ray  angular  flux 
spectra  calculated  with  ANISN  at  17 
angles  for  a  12.2-15  MeV  neutron 
sources  at  a  range  of  2400  meters. 
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Ratio  functions  as  defined  earlier  were  formed  into  a  data  base 

I 

R(Eg,  E,  r,  m)  for  modeling  and  are  illustrated  in  the  representative 
plots  of  Figs.  36  and  37  for  neutrons  and  secondary  gamma-rays, 
respectively,  at  300  meters  from  a  14  MeV  source.  These  curves  show 
clearly  the  smoothing  effect  as  a  function  of  energy  inherent  in  the  ratio 
function.  For  neutrons,  R  is  unity  to  a  good  approximation  for  energies 
below  0. 1  MeV.  Above  1  MeV,  considerable  structure  was  present  in 

I 

R  which  made  it  difficult  to  fit  well  with  simple  functions.  The  slowly 

i  1 

varying  nature  of  the  angular  flux  vs  distance  (Fig.  35)  indicated  that  the 

I 

ratio  could  be  parameterized  more  efficiently  by  using  distance  as  an 
independent  variable.  Thus,  the  following  representation  was  used  to 
model  the  ratio  function  of  neutrons  for  each  source  energy: 


R(E,  r,n)  =  1.0  for  E<0.1  MeV 


RE(r,  fi)  =  exp  Aj(fi)  +  A2(/i)r+A3(M)r' 


for  E>0. 1  MeV 


where 

A(m)  =  B^expfBg#*)  for  -1.  0  <  n  <  0.  6178 

=  exp  (C  j+Cg  ^  *or  0.  6178  <  f*<  1.0 

t 

This  procedure  increased  the  total  number  of  neutron  coefficients  by  less 
than  a  factor  of  2  to  represent  17  times  the  amount  of  data.  Thus,  the 
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Fig.  36.  Neutron  ratio  function  for  14  MeV  neutron  source  at 
300  meters. 
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entire  energy-angle  data  base  for  neutrons,  representing  over  10  flux 
values,  required  about  2000  coefficients. 

The  ratio  function  for  secondary  gamma-rays  (Fig.  37)  varied 
more  smoothly  with  energy  and  was  fit  versus  energy  for  each  source  by 
a  five -order  polynomial  whose  coefficients  themselves  were  sub-poly¬ 
nomials  in  distance  and  angle.  These  parametric  fits  essentially  doubled 
the  number  of  coefficients  in  the  secondary  gamma-ray  model  for  a  17- 
fold  increase  in  flux  information. 

The  accuracy  of  the  neutron  and  secondary  gamma-ray  angular 
flux  model  is  difficult  to  estimate  because  of  the  inherent  limitation  in  the 
raw  data  base.  Differences  between  the  parametric  model  and  the  raw 
ANISN  data  base  can  be  large  for  the  small  fluxes  at  high  energies  and  in 
the  backward  directions,  however,  the  raw  data  in  this  region  are  probably 
uncertain  by  a  factor  of  2.  The  agreement  between  the  model  and  the 
"smoothed"  angular  flux  data  base  is  of  course  much  better.  The  most 
meaningful  estimate  of  the  adequacy  of  the  angular  flux  model  is  probably 
obtained  by  comparing  partially  integrated  quantities  with  the  raw  ANISN 
data  base.  Thus,  the  total  fluences  integrated  over  energy  and  angle 
generally  agreed  within  5%;  fluence  distributions  integrated  over  energy 
aid  o  agreed  to  about  10%;  and  fluence  distributions  integrated  over  angle 
agreed  to  about  20%.  Differential  angular  flux  data  calculated  by  the 
parametric  model  are  shown  in  Figs.  38  and  39,  for  neutrons  and  secondary 
gamma-rays  at  600  meters  from  a  14  MeV  source.  These  curves 
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Fig.  39.  Secondary  gamma-ray  energy -angle 
distribution  calculated  by  ATR  at  600 
meters  from  a  14  MeV  neutron  source. 
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demonstrate  the  desired  structure  and  the  systematic  behavior  charac¬ 
teristic  of  deep-penetration  transport. 

The  extension  of  the  angular-flux  modeling  to  prompt  photon 
transport  using  a  ratio  function  is  expected  to  be  a  less  tedius  job  due 
to  the  lack  of  energy  structure  in  the  scalar  flux  shapes.  Preliminary 
estimates  indicate  that  the  angular  dependence  of  prompt  photon  trans¬ 
port  will  require  approximately  1000  coefficients  each  for  the  low  and 
high-energy  ranges. 
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IV.  THE  ATR  CODE 

The  parametric  models  described  in  the  previous  section  have 
been  used  as  the  basis  for  constructing  ATR  -  a  computer  code  for 
efficiently  calculating  the  air  transport  of  radiation  for  a  variety  of  differ¬ 
ent  problems  at  all  altitudes  in  the  atmosphere.  Mass  scaling  is  used  to 
transform  the  parameterized  results  for  homogeneous  air  from  an  air 
density  of  1.11  mg/cm^  to  other  air  densities.  Corrections  are  also 
included  for  air-ground  interface  and  exponential  air  density  effects.  The 
code  has  been  structured  to  allow  the  user  a  maximum  degree  of  flexibility 
in  calculating  different  problems,  with  a  minimum  number  of  restrictions 
on  input  data  and  format.  The  presently  completed  version  of  ATR  calculates 
problems  for  neutron  and  secondary  gamma-ray  transport.  Prompt 
photon  transport  for  high  and  low  energy  photons  will  be  included  in  the 
near  future  when  modeling  is  completed  on  photon  angular  distributions. 

4. 1  Problem  Geometry  and  Scaling  Laws 

The  problem  geometry  for  ATR  is  illustrated  in  Fig.  40, 
where  the  geometry  coordinates  are: 

H  =  source  altitude 
s 

H^,  =  target  altitude 

R  =  slant  range 
s 

RH  =  horizontal  range 
9  -  slant  angle 
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The  fundamental  calculation  in  ATR  is  to  determine  tin*  free- 


field  distribution  lor  a  givtm  problem  geometry  and  source  spectrum. 

In  order  to  calculate  such  problems  for  the  actual  atmosphere,  two 
assumptions  are  made. 

1.  The  transport  of  radiation  between  source  and 
target  can  be  approximated  by  homogeneous  air  trans¬ 
port  results  for  a  density  equal  to  the  average  density 
between  the  source  and  target  altitudes. 

2.  The  corrections  due  to  air-ground  interface 
and  exponential  air  density  effects  can  be  included  as 
perturbations  to  homogeneous  air  results. 

The  general  cal»~alational  scheme  in  ATR  is  therefore  to  scale  transport 
results  from  the  parametric  model  density  of  1.11  mg/cm  to  the  average 
density  of  the  problem  and  then  to  include  correction  factors  for  the 
effects  of  an  air-ground  interface  and  exponential  air  density.  These 
correction  factors  were  essentially  negligible  for  altitudes  from  1-20  KM. 
Density  scaling  laws  for  the  Boltzmann  equations  are  reviewed  in 
Appendix  B  and  show  that  the  angular  flux  0t(E,  il,r)  in  the  transformed 
system  of  density  Pt(rt)  is  related  to  the  flux  in  the  original  system  at 
density  p(r)  by 

0(E,fl,r)  =K20t(E,fl,  rt) 

where 
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Kr  and  Pt(rt)  =  p(r)/K 


i. e. ,  for 


ptrt 


P  r 


2 

This  is  the  well  known  result  that  the  Anr  flux  depends  only  on  the 

mass  of  air  through  which  the  neutrons  and  photons  are  being  transported. 

2 

The  ATR  parametric  model  is  based  on  An r  quantities  as  a  function  of 

2 

penetration  distance  in  grams/cm  .  The  total  mass  of  air  between 

source  and  target  for  a  given  problem  is  calculated  from  R  !>  and  used 

s 

to  obtain  4/rR^  quantities  from  the  parametric  data  base,  i.e. , 
s 


For 


4  nR  <P  =  A  nR  $ 

S  S  0  0 


R  P  =  R  p 

S  0  0 


where 

P  -  average  air  density  of  the  problem 

o 

Pq  air  density  of  ATR  data  base  (1.11  mg/cm  ) 

<P  radiation  free-field  calculated  by  parametric 
equations  at  density  PQ  and  distance  Rq 

<P  -  radiation  free-field  of  the  problem  geometry 
s 

The  average  air  density  of  the  problem  is  calculated  from 
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where  p(Z)  is  the  variation  in  the  atmospheric  density  versus  altitude 
above  sea  level. 

The  atmospheric-density  was  calculated  for  altitudes  from 

(74) 

0-300  KM  based  on  a  model  developed  by  Dupree.  Dupree’s  model 

consisted  of  a  product  of  32  exponential  terms  and  was  used  for  altitudes 

from  0-100  KM.  Above,  100  KM,  an  empirical  curve  fit  was  used.  The 

entire  density  model  was  based  on  data  listed  in  the  U.  S.  Standard 
(78) 

Atmosphere.'  Equations  used  for  the  model  were: 


-Pj(h-hM) 

P(h)  =  * - 

?o 


-pi(hi-hi.1)i 


0<h<100  KM 


for  h.  .<h<h. ,  with  h  =0  and  3<H32.  The  product  is 
I-l_  i  o 

deleted  for  1=2,  and  PQ  is  the  sea  level  density  value. 
p(h)  =  exp  [Cj+Cgh+Cgh2]  100  KM<  h  <  130  KM 

r  2^t 

p(h)  =  exp  I  Cj+Cph+CJi  + -  130  KM<  h 

1  L  6  h-99J 
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4.2  Air-Ground  Interface;  Corrections 


Corrections  for  air-ground  interface  effects  were  based  on  the 

(75) 

first-last  collision  model  of  French.  In  this  model  factors  F(Hg) 
and  F(H,p)  are  derived  for  correcting  infinite  air  results  according  to 

0(Hs,Ht,  R)  f(hs)f(ht)$o(r) 

where  d>  (R)  is  the  infinite  air  flux, 
o 

The  model  estimates  F(Hy)  by  calculating  the  effective  fraction  of  first 
collisions  about  the  source  as  a  function  of  source  altitude.  F(H,p)  is 
determined  by  calculating  the  effective  number  of  last-collisions  in  the 
vicinity  of  the  target.  The  model  was  developed  for  fast  neutrons  and 
treats  neutron  reflection  from  the  ground  by  an  albedo  approach.  The 
correction  factors  are  listed  in  Table  II  as  a  function  of  MFP  from 
the  ground. 

Comparisons  of  the  model  prediction  with  Straker's  two-dimen- 

(12) 

sional  transport  calculations  of  dose  vs.  distance  showed  that  the  model 
was  more  accurate  in  correcting  from  one  source  height  to  another  than 
correcting  infinite-air  results.  This  was  also  true  for  correcting  sec¬ 
ondary  gamma-rays  where  the  model  predictions  alone  gave  results 
differing  from  Straker’s  calculations  by  as  much  as  a  factor  of  five. 
Correcting  air-ground  data  from  one  source  altitude  to  another,  however, 
gave  agreement  to  about  25%. 


66. 


TABLE  n. 

Air-ground  correction  factors  from  the 
first-last  collision  model  as  a  function 
of  MFP  from  the  ground.  > 


H(MFP) 

F(HS) 

F(Hr) 

0 

0.434 

0.01 

0.466 

0.717 

0.02 

0.478 

0.728 

0.05 

0.506 

0.750 

0.10 

0.546 

0.781 

0.25 

0.632 

0.841 

0.50 

0.724 

0.902 

0.75 

0.798 

0.934 

1.00 

0.862 

0.953 

1.50 

0.941 

0.  977 

2.00 

0.972 

0.987 

2.50 

0.987 

0.990 

0.993 

0.995 

3.00 

0.993 

3.50 

0.996 

4.00 

0.997 

0.997 

H>4. 0 

1.00 

1.00 

On  the  basis  of  the  above  considerations,  a  scheme  was  developed 
for  obtaining  air-ground  corrections  to  infinite  air  results  by  using  the 
first-last  collision  model  to  extrapolate  Straker’s  dose  data  to  other 

I  i 

altitudes,  i.e.,  for  H,p<100  meters: 


C(Hg,  HT,Rg) 


F(Hs)F(Ht) 

F(Hs)F(Ht) 


DInf.  Air(HS’  HT’  RSI 


where 


Ilg  -  15  meters 
H^,  =  1  meter  ' 

D  -  Henderson  dose  vs  distance  calculated  by  Straker 
0  at  Hg  and  H^,  , 

Hg,  Ht  and  Rg  are  the  problem  coordinates 

3 

PQ  =  1.11  mg/cm 

I 

p  average  air  density  in  problem  geometry 


The  above  correction  factor  is  used  in  ATR  to  correct  all  differential 
results  obtained  by  using  mass  scaling  of  the  parametric  model  for 
target  altitudes  below  100  m.  The  two  main  assumptions  in  the  above 
correction  scheme  are: 

1.  The  spatial,  energy  and  angular  distributions 
of  radiation  near  an  air-ground  interface  can  be  approxi¬ 
mated  from  infinite  air  results  with  a  shift  in  intensity  only. 


2.  The  intensity  shift  is  proportional  to  ratio 
of  the  air-ground  dose  to  the  infinite  air  dose  cal¬ 
culated  for  transport  through  the  same  mass  of  air. 

For  target  altitudes  above  100  meters,  the  infinite  air  results  are 
corrected  using  the  first-last  collision  model  alone,  i.e.,  for  HT>100 
meters: 

c(hs,ht,rs)  =  F(Hg)F(HT)  . 

The  effect  of  the  air-ground  corrections  is  illustrated  in  Fig.  41,  which 

shows  dose  vs  distance  curves  for  the  Hiroshima  device  obtained  using 

(77) 

ATR  and  compared  with  the  empirical  results  of  Auxier  et.  al.  The 
ATR  results  with  air-ground  corrections  are  approximately  50%  lower 
than  scaled  infinite  air  results,  and  are  in  fair  agreement  with  Auxier's 
values. 

4.3  Exr*'-  ntial  Air-Density  Corrections 

Corrections  for  non-uniform  air  density  effects  were  obtained 
by  extending  the  first-last  collision  model  to  include  the  density  varia¬ 
tion  versus  altitude  in  the  numerical  integration  scheme.  The  first- 
collision  correction  factor  for  source  altitude  was  obtained  by  numerical 
integration  of  the  equation  below.  ' 

F(Hg) .  £  f  r  wmL y(z)^xdz 

a  *  J  J  x2+z2 

0  -HS 
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Fig.  41.  Neutron  dose  from  Hiroshima  device  (12.5  KT 
yield,  579  m  height-of-burst). 
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where  Mq,  is  the  total  cross  section  in  air  in  cm  /gm.  W(Z)  is  a 
weighting  factor  given  by  F(Hg+Z)  to  account  for  multiple  scattering 
and  is  obtained  by  an  iteration  of  the  above  equation. 

The  additional  term  used  to  account  for  reflection  in  the  air- 
ground  case  is  zero  here  since  neutrons  are  not  reflected  from  the  top 
of  the  atmosphere. 

The  last-collision  correction  for  target  altitude  was  obtained 
from  numerical  evaluation  of  the  above  integral  without  the  weighting 
function.  The  numerical  results  for  F(Hg)  and  F(HT)  are  listed  in 
Table  III  as  a  function  of  altitude  in  terms  of  mean  free  path  to  the 
top  of  the  atmosphere.  These  correction  factors  were  applied  uni¬ 
formly  to  the  scaled  infinite-air  results  in  ATR  for  both  neutron  and 
secondary  gamma-rays.  The  validity  of  the  model  was  checked  against 
the  two-dimensional  Monte  Carlo  calculations  of  Keith  at  the  two 
source  altituder  reported  of  85  KFT  and  110  KFT.  The  first-last 
collision  model  reduced  the  infinite-air  flux  by  factors  which  approxi¬ 
mated  the  average  behavior  of  the  Monte  Carlo  results,  but  did  not 
predict  the  same  relaxation  lengths.  At  large  distances  where  the 
model  is  expected  to  be  a  good  approximation,  the  correction  factors 
were  in  fairly  good  agreement  with  Monte  Carlo  results.  At  a  co- 
altitude  distance  of  50  gm/cm  ,  for  example,  the  correction  factors 
were  .  79  and  .  23  respectively,  for  a  14  MeV  source  at  altitudes  of  85 
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TABLE  III 


Exponential  air  correction  factor  from  the  first-last 
collision  model  as  a  function  of  MFP 
from  the  top  of  the  atmosphere. 


H(MFP)  F(HS)  F(HT) 


0 

0 

0 

0.010 

0.067 

0.180 

0.030 

u92 

0.229 

0.060 

0.133 

0.300 

0.090 

0.166 

0.350 

0.140 

0.210 

0.408 

0.210 

0.271 

0.482 

0.320 

0.373 

0.585 

0.490 

0.491 

0.676 

0.770 

0.642 

0.793 

1.200 

0.808 

0.890 

2.000 

0.920 

0.967 

3.100 

0.975 

0.980 

4.900 

0.985 

0.990 

7.900 

0.996 

0.998 

H>7.  9 

1.000 

1.000 

KFT  and  110  KFT;  these  compared  well  with  the  corresponding  ratios  between 

the  non-uniform  and  uniform  air  results  determined  from  Keith's  curves  and 

were  0. 69  and  0. 23.  Thus,  the  first-last  collision  model  seems  to  predict 

the  average  effect  of  non-uniform  air  on  neutron  transport.  Secondary  gamma- 

ray  data  were  not  available  for  comparison. 

The  effect  of  the  non-uniform  air  correction  in  ATR  is  demonstrated 

by  the  isodose  curves  in  Fig.  42  which  were  calculated  at  a  source  altitude 

of  90  KFT  for  an  isodose  value  of  2500  rads.  The  source  was  the  unclassified 

26 

thermonuclear  spectrum  of  ORNL  4464  with  a  yield  of  10  neutron.  The  non- 
uniform  air  correction  result  in  slant  ranges  which  are  smaller  than  for 
scaled  homogeneous  air  since  the  corrected  flux  curve  is  shifted  to  smaller 
flux  values. 

4.4  ATR  Control  Commands 

The  list  of  control  commands  for  setting  up  an  ATR  run  J  j  given  in 
Fig.  43.  For  the  most  part,  the  order  of  control  commands  is  immaterial 
for  any  given  set  of  problems,  i.  e. ,  between  *EXC  commands.  The  initiali¬ 
zation  routine  contains  some  calculations  on  parameters  which  can  be  evaluated 
on  the  basis  of  other  parameters  as  well  as  some  storage  initialization.  When 
a  command  is  read  in  and  interpreted,  the  corresponding  subroutine  is  called 
to  interpret  and  properly  place  all  of  the  parameters  following  the  command. 
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Fig.  42.  Isodose  curves  calculated  with  and  without 
non-uniform  air  correction  for  a  source 
altitude  of  90  KFT. 
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Errors  in  control  command  format  or  structure  are  noted  immediately  but 
no  other  visible  action  is  taken  by  the  code  until  the  *EXC  command  at  which 
point  output  is  generated  which  may  be  just  a  few  numbers  or  hundreds  of 
pages,  depending  on  the  input  parameters. 

The  code  also  contains  a  set  of  subroutines  which  permit  an  essentially 
format-free  placement  of  numbers.  The  reading  of  values  is  upward  com¬ 
patible  in  that  real  values  need  not  contain  decimal  points.  Standard  FORTRAN- 
type  real  numbers  in  scientific  notation  are  acceptable  with  or  without  the 
E;  i.  e. ,  the  number  400  can  be  written  as  400  or  400.  or  4.  E+2  or  4.  E2 
or  4. +2.  The  blank  character  is  used  as  a  number  delimiter,  thus  numbers 
must  not  contain  blanks.  Any  number  of  blanks  may  separate  numbers  and 
a  series  of  values  may  be  continued  on  subsequent  cards,  but  a  number  can¬ 
not  be  split  between  two  cards.  The  major  restriction  is  that  since  the  real 
number  and  fractional  parts  as  integers,  neither  the  whole  number  nor  the 
fractional  part  of  any  real  number  can  exceed  the  integer  capacity  of  the 
particular  computer. 

Control  commands  and  their  associated  parameter  values  are  cumu¬ 
lative  so  that  even  after  the  results  of  a  particular  problem  have  been  dis¬ 
played,  the  parameter  values  and  flags  remain  the  same  for  the  next  problem 
unless  specified  anew  with  another  command  of  the  same  type.  Consequently, 


76. 


from  one  problem  to  the  next,  only  the  most  essential  parameters  need  to 
be  specified.  This  feature  allows  multiple  runs  to  be  made  easily  which 
require  only  a  few  parameter  changes  to  be  made  for  each  problem.  Thus, 
the  following  sequence  of  commands  is  allowed: 


*EXC 

*EXC,  4PIRSQ 

2 

This  sequence  executes  the  particular  problem  first  without  4»rR  and  then 
2 

with  4/rR  included  as  a  factor  in  the  output  quantities.  The  only  command 
exempt  from  this  is  the  *CONSTRAINT.  In  order  to  clear  out  unwanted 
parameters  the  *STOP  command  is  used  after  which  the  problem  must  be 
completely  respecified. 

4. 5  Input 

The  list  of  input  control  commands  which  were  listed  in  Fig.  44 
will  be  explained  in  detail  in  this  section. 

4.  5.  i  *Z-SOURCE(I) 

Z  =  N,  G,  X  and  denotes  neutron,  prompt  gamma  or  X-ray 
source,  respectively,  and  has  the  same  meaning  for  other  source- 
related  commands.  I  =  1 ,  2,  3,  4,  5,  or  6  where  1-3  stands  for  internal 
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1.  *Z-S0URCE(1) 

2.  *Z-SVAL,  units,  values 

3.  *Z-EVAL,  units,  values 

4.  *Z-NORM  value 

5.  +  Z-YIELD  value 

0.  +XX,  units,  value(s) 

7.  *E-REGRP/Z/,  units,  values 
3.  *A-REGRP/Z/,  units,  values 

9.  *  PRINT/ Z/ (II  12  .  .  .  ) 

10.  *WRITE/Z/ (II  12  .  .  .  ) 

11.  *C0NSTRAINT/Z/ (XX  I  Value),  units 

12.  *CXX,  units,  value(s) 

13.  *FLUXWT/Z/,  units,  values 

14.  *TITLE  n 

15.  +DOSE/Z/ 

16.  +EXC,  4PIRSQ 

17.  +STOP 

18.  +FIN 

19.  ♦GROUND,  units,  value 


Fi|/.  44.  ATR  input  commands. 
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source  selections,  4  denotes  external  histogram  values,  5  denotes  exter¬ 
nal  point  values  and  6  specifies  fission/fusion  source.  The  internal  source 
options  taken  from  ORNL  4464  are  fission  spectrum  (1=1),  thermonuclear 
(1=2)  and  12-15.2  MeV  monoenergetic  band  (1=3)  for  neutrons.  Options  4 
and  5  require  the  use  of  the  *Z-SVAL  command  to  read  in  source  values 
and  a  *Z-EVAL  to  read  in  energy  values.  If  the  internal  source  energy 
structure  is  used  then  only  the  *Z-SVAL  command  is  necessary.  The  Value 
is  the  fraction  (between  0  and  1)  of  12-15.  2  MeV  neutrons  in  the  internal 
fission  spectrum  and  it  is  required  for  option  6  only. 

Example:  *N-SOURCE(4) 

4.5.2  *Z-SVAL,  Units*,  Values** 

Units  =  PER  MEV,  PER  KEV  or  PER  GROUP  and  defines  the  source 
units.  If  units  are  not  specified,  the  default  choice  is  set  internally  to 
PER  GROUP.  Values  =  source  values  to  be  read  in.  These  values  must 
correspond  to  the  energy  values  read  in  by  the  *Z-EVAL  command.  The 
order  is  assumed  to  be  from  low  energy  to  high  energy. 

Example:  *N-SVAL  .01  10  33  100  20. 

+If  units  are  designated  with  any  of  the  commands,  they  must  be  delimited 
by  commas;  otherwise,  default  unit  designations  are  used  and  commas 
must  not  appear. 

++Va!ues  must  be  separated  by  at  least  one  blank  and  may  be  continued  on 
the  next  card. 
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4.5.3  *Z-EVAL,  Units,  Values 


Units  =  MEV  or  KEV.  A  default  in  units  specifies  MEV  for  N  and 
G,  and  KEV  for  X.  Values  =  source  energy  boundary  values  if  (4)  is  spec¬ 
ified  in  4.  5. 1,  or  source  energy  point  values  if  (5)  is  specified  in  4.  5. 1. 

The  order  is  from  low  to  high  energy  values. 

Example:  *N-EVAL,  KEV,  3.35  12  50.5  120  170  235 

4.5.4  *Z-NORM  Value 

Value  =  source  normalization  value  in  units  of  particle/KT.  If 
this  option  is  not  specified  then  the  source  spectrum  will  not  be  normalized. 
Example:  *N-NORM  2. 

4.5.5  *Z -YIELD  Value 

Value  =  source  yield  in  KT.  The  total  source  output  is  the  product 
of  the  YIELD  value  and  the  NORM  value.  Default  value  is  1.0. 

Example:  *N- YIELD  1 . 5E+23 

4.5.6  *XX,  Units,  Value(s) 

XX  =  one  of  RH,  RS,  HT,  HS  and  AN.  It  denotes  the  horizontal 
range,  slant  range,  target  height,  source  height  and  slant  angle  respec¬ 
tively.  Units  -  one  of  M,  KM,  MILE,  YD,  KFT,  and  FT  for  the  distance 
and  one  of  DEG,  RAD,  COS  for  the  angle  specification.  Distance  default 
unit  is  M  and  angle  default  unit  is  DEG.  Values  =  up  to  50  values  of  the 
geometry  parameter.  These  values  may  be  read  in  one  at  a  time 
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separated  by  blanks  or  in  the  format  n^(n)ng  signifying  values  from  n^  to 
n g  in  steps  of  n.  In  this  case  of  the  latter  must  be  arithmetically 
greater  than  n^ .  Three  geometry  values  must  be  specified  for  a  meaning¬ 
ful  geometry  configuration. 

Two  of  the  possible  geometry  configuration  definitions  result  in 
ambiguities.  One  is  when  RH,  RS  and  HT  are  specified.  There  is  no 
inherent  knowledge  in  this  specification  of  whether  HS  should  be  placed 
above  or  below  HT.  In  order  to  resolve  the  ambiguity  the  characters 
"+"  or  should  be  used  with  the  *HT  command  to  indicate  that  HS  is 
to  be  placed  above  or  below  HT,  respectively.  The  other  ambiguous 
configuration  occurs  when  RH,  RS  and  HS  are  specified.  In  this  case 
the  4  or  -  should  be  used  with  the  *HS  command  to  place  HT  above  or 
below  HS  respectively.  If  none  of  these  characters  are  specified  then 
+  is  assumed. 

There  is  a  special  units  option  (GM)  on  RS  which  is  allowed  if  the 

following  geometry  configuration  is  specified:  HS,  HT,  RS.  In  that  case 

the  values  following  the  "*RS,GM,"  command  are  interpreted  as  units 
o 

of  gm/ cm  of  the  slant  range.  The  GM  units  option  is  restricted  to  this 
configuration  only. 

Examples:  *HT+,  KM,  1. 

*RS,  KFT,  .5(.  5)10 
*RH,  KM,  10 
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4.5.7  *E-REGRP/Z/,  Units,  Values 

Z  =  N,  G,  X,  NG,  GGfor  neutron,  prompt  gamma,  X-ray,  sec¬ 
ondary  gamma  arid  combined  prompt  and  secondary  gammas,  respectively. 
Units  have  the  same  options  and  meanings  as  in  4.  5.  3.  Values  =  energy 
boundary  values  (low  to  high)  which  will  be  used  to  regroup  output  results. 
Example:  *E-REGRP/NG/  .  1  .5  1.  2.  3.  5.  10 

4.5.8  *A-REGRP/Z/,  Units,  Values 

Z  has  the  same  options  and  meanings  as  in  4.  5.  7.  Units  has  the 
same  meaning  as  the  angle  specification  in  4.  5.  6.  Values  -  angle  bound¬ 
ary  values  (low  to  high  cosine  values)  which  will  be  used  to  regroup  output 
results. 

Note  that  if  angles  are  specified  in  degrees  for  this  command  then 
the  angle  values  must  be  entered  high  to  low. 

Example:  *A-REGRP/N/,  COS,  -1.  -.5  0  .25  .5  .75  1.0 

4.5.9  * PRINT/Z/ (1 1  12  - --) 

Z  has  the  same  meanings  and  options  as  in  4.5.7.  II,  12,  etc.  - 
integer  values  from  0  to  16,  and  will  yield  the  following  results  on  the 
print  output  unit: 

0  s  summary 

1  -  number  fluence 

2  number  current 

3  energy  fluence 


82. 


4  =  energy  current 

5-12  -  internal  dose  responses 

13  ^  external  flux  weight  as  read  in  by  *FLUXWT  command 

14  not  used 

15  =  results  from  *DOSE  command 

16  =  results  from  "CONSTRAINT  commands 

If  any  of  these  appear  with  a  negative  sign,  then  the  angle-energy  distri¬ 
bution  output  will  be  suppressed  for  the  particular  case  (applied  to  1-13). 
The  internal  dose  responses  taken  from  ORNL-4464  correspond  to  the 
following  specifications: 

Neutron^ 

5  -  Henderson  Tissue  Dose 

6  -  Snyder-Neufeld  Dose 

7  -  Tissue  Kerma 

8  -  Mid- Phantom  Dose 

9  -  Concrete  Kerma 

10  -  Air  Kerma 

11  -  1  MeV  Silicon  Equivalent 

12  -  Ionizing  Silicon  Kerma 

Example:  *PRINT/NG/  (0  1  -3  2  6  -13) 

Options  1-4  are  defined  in  Section  4.  6.  The  summary  (option  0) 
includes  total  integrated  quantities  of  currents,  fluences,  external  flux 
weight,  average  energy,  Henderson  tissue  dose  and  ionizing  silicon  dose. 


Gamma* 

5  -  Henderson  Tissue  Dose 

6  -  Concrete  Kerma 

7  -  Air  Kerma 

8  -  Silicon  Kerma 

9-12  -  not  used 


^All  responses  are  in  units  of  rad/(particle/cm^). 
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Options  15  and  16  are  automatically  specified  by  the  corresponding  com¬ 
mands  and  need  not  be  explicitly  specified  by  the  *PRINT  command. 
4.5.10  *WRITE/Z/(I1  12  — ) 

All  of  the  options  are  the  same  as  in  4.  5.  9. 

Options  15  and  16  are  automatically  defined  for  *  PRINT  but  must 
be  explicitly  specified  for  *WRITE. 

Unit  numbers  are  set  up  in  a  DATA  statement  where  they  can  be 
easily  modified  for  any  particular  computer  system.  The  current  values 
are: 

5:  input  unit 
6:print  unit 
8: write  unit 

The  *WRITE  command  is  included  so  that  production  runs  generated  by 
ATR  can  be  saved  on  magnetic  tape,  drum  or  disk  for  future  use. 

Example:  *WRITE/N/  (2  6  7  -8  15  16) 

4.  5.  H  »CONSTRAINT/Z/(XX  I  Value),  Units 

Z  has  the  same  meanings  and  options  as  in  4.  5.7  and  I  has  the 
same  meanings  and  options  as  in  4.  5.  9,  but  only  options  1-13  are  mean¬ 
ingful.  Value  is  the  fluence,  current  or  dose  value  for  which  the  con¬ 
straint  geometry  is  to  be  evaluated.  This  command  only  makes  sense 
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for  response  values  that  are  not  multiplied  by  4wR  .  Units  has  the  same 
meanings  and  options  as  in  4. 5.  6  and  it  is  a  specification  defining  the 
geometry  units  printed  in  the  output. 

The  ^CONSTRAINT  command  calculates  the  value  of  coordinate 
XX  which  corresponds  to  a  predetermined  response  level  for  given  values 
of  two  other  geometry  coordinates  specified.  Version  1  of  the  code  is 
restricted  to  XX  =  RS. 

Example:  *CONSTRAINT/N/  (RS  8  5E-20),  MILE 

4.5.12  *CXX,  Units,  Value (s) 

*CXX  defines  the  fixed  coordinates  and  their  values  for  a  con¬ 
straint  calculation.  All  parameters  have  the  same  meaning  and  options 
as  in  4.  5. 6. 

Two  *CXX  commands  combined  with  a  ^CONSTRAINT  command 
define  a  constraint  problem;  of  course,  care  must  be  taken  to  specify 
consistent  and  valid  geometry  configurations. 

Examples:  *CHS  1000 

♦CAN,  DEG,  -90  (10)  90 

4. 5. 13  *FLUXWT/Z/,  Units,  Values 

Z  has  the  same  meanings  and  options  as  in  4.  5.7.  Units  = 
arbitrary  (up  to  16  characters)  unit  specification  that  will  appear  on 
the  summary  output.  Values  =  flux  weights  to  be  used  for  weighting 
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I  ho  ATR  energy  spectrum  and  should  orrespond  to  the  internal  energy 
structure  (low  to  high). 

Examples: 

♦FLUXWT/NG/,  CM**2  ERGS/ GRAM,  5.9-8  1.2-8 
8.-9  1.2-8  1.7-8  2.4-8  3.3-8  4.1-8  5.-8  6.15-8 
7.15-8  8.4-8  9.8-8  1.2-7  1.5-7  1.8-7  2.2-7  2.65-7 

4.5.14  ♦TITLE  n 

n  may  be  up  to  74  characters  and  is  used  as  a  problem  title  to 
identify  the  output. 

Example:  *TITLE  AIR  TRANSPORT  PROBLEM 

4.5.15  *DOSE/Z/ 

This  command  automatically  prints  dose  values  vs  the  running 
coordinate  value  specified.  Z  has  the  same  meanings  and  options  as  in 
4.5.7. 

Examples:  ♦DOSE/NG/ 

♦DOSE/N/ 

4.  5. 16  *EXC ,  4PIRSQ 

This  command  is  used  to  execute  the  problem  defined  by  prior 

commands.  If  the  comma  appears  after  *EXC  then  resulting  values  will 

2 

be  multiplied  by  4?rR  where  R  is  the  slant  range. 
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4.5.17  *STOP 


This  command  clears  the  internal  flags  so  that  the  next  problem 
can  be  specified.  All  previous  specifications  are  eras'd  by  this  command, 
and  it  would  normally  follow  a  *EXC. 

4.5.18  *FIN 

This  command  terminates  the  program. 

4.5.19  *  GROUND,  Units,  Value 

Units  has  the  same  meanings  and  options  as  in  4.  5. 6.  Value  is 
the  ground  level  relative  to  sea  level  for  the  ATR  problem  geometry.  If 
this  command  is  not  specified  then  sea  level  will  be  used  for  the  ground 
level. 

4.  6  ATR  OUTPUT  AND  SAMPLE  RUNS 

Five  basic  quantities  can  be  calculated  from  ATR  for  each 
source-particle  type: 

2 

Number  Fluence:  <p  (r,  E,n)  particles/MeV  •  cm  •  steradian 

2 

Number  Current:  pw(r,  e,p)  particles/MeV  •  cm  •  steradian 

2 

Energy  Fluence:  E<;i(r,  E,m)  particle-MeV/MeV-  cm  •  steradian 

2 

Energy  Current:  Efi<ff(r,  E,p)  particle-MeV/MeV- cm  -steradian 

DOSE:  D(E)  p(r,  E,/x)  rads/MeV*  steradian 


i 
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The  ATR  control  commands  determine  the  type  of  output  desired  by  the 
user,  can  consist  of  the  differential  quantities  above,  or  any  combination 
of  the  above  integrated  over  energy  or  angle.  A  series  of  three  sample 
ATR  runs  was  made  on  a  Univac  1108  computer  to  illustrate  use  of  most 
of  the  input  commands.  The  computer  print-out  is  reproduced  in  this 
section  and  illustrates  some  of  the  output  quantities  available  from  ATR 
and  the  corresponding  input  required  for  each  run. 
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SAMPLE  PROBLEM  1 
Page  1 


VISION  1,  Of  ATRNXC 

♦  HUE  SAMPLE  ATR  PROHLEm  i 
K-VjUKCU?) 

•N-NLW-  1,0*2? 

♦  N-VIE1.0  100 

i  1 

•■‘bil'M 

♦Hi .-1,0.621 

♦  I,  .5  >  8  f  l  I  o ,  5  6 

♦  H.‘»*T/n/.R»OS#  I  1  l  1  I  1  1  1  1  1  1  1  1  1  I  1  1 
•PKINt/N/IO) 

♦  PH  TO  /MU/  (01 
»C  XC. 
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AIR  PROHLEM  RUMMER 


SAMPLE  AIR  PROHLEm  1 


ffGTRON  SOURCE  INTERNAL  THERMONUCLEAR 


NORMALIZATION*  1.000*2?  NEUTRON  /KT, 
TOTAL  OUTPUT*  1.000*29  NEUTRON 


SOURCE  SPCCTRUM 


ENERGY (Mf  v  ) 

N/KT  N/K  T  *MEV 

1,07-05- 

2.90-05 

0.00 

0,00 

?. 90-05- 

1,01-09 

?.00*?1 

2.78*25 

1. Ol-oa- 

5. OS-flU 

2.90*2? 

0,98*25 

5.8J-00- 

3.T5-0S 

l.??*?l 

9,91*25 

3.35-03- 

1,11-01 

3.65*23 

3.39*29 

1.11-01- 

5.S0-01 

1 .02*23 

2.32*23 

5. 50-01- 

1. 11*00 

8.50*?? 

1.52*23 

1,11*00- 

1.81*00 

6, 20*22 

8.81*22 

l ,*5*00- 

2.35*00 

2. 80*22 

5.38*22 

TtFLB*  1. 000*0?  KT 


ELO*  l.o 

ENERGYEmev  > 

N/KT  1 

m/kT*me  V 

2,35*00- 

2,98*00 

5.00*21 

9.55*22 

2.98*00- 

3,01*00 

1.90*22 

3.95*22 

3,01*00- 

9,07*00 

2.80*22 

2.95*22 

9.07*00- 

9,97*00 

1.70*22 

l.M9*2? 

9  c 97*00- 

8.18*00 

1.80*22 

1.29*22 

6.36*00- 

8,19.00 

1.97*22 

8,03.21 

8,19*00- 

1,00*01 

T .91*22 

7.79*21 

1.00*01- 

1,22*01 

2.58*22 

1.18*22 

1.22*01- 

1.50*01 

7.08*22 

2.52*2? 

ATM  PROHLfM  NUMHFR  1  SAMPLF  ATR  PROBLEM  1  •  •  *  * 


GROUND  LEVEL 
•MCR17.  RANGE 
SLANT  RANGE 
TARGET  ALT, 
SOliRCt  ALT, 
•SLANT  ANGLE 


HH* 

RS« 

HT» 

HS* 

AN* 


,305«M, 

1,999«M, 

1.999KM, 

,999«m, 

l.oqoKM, 


36.791GM/CM**?, 
2 1  5 «  79 1  G*VCM**2 , 
?15.79«G“/Cm*»?. 
1 1  5.?73Gm/Cm**2, 
11  J.1J7Gm/Cm*»?, 


•,01 70EGREES  (COS*  1,00000) 


l .OOOKE  T, 
6.560KET. 
6.560*8  T  , 
3.279NNT, 
J,?81KFT, 


•CALCULATED  FROM  OTHER  COOROlNAlES 


.190VJLFS 

t.?'<?»rLES 

l.?R?R!LES 

,6?1"ILES 

,fr?l“ILFS 


NEUTRON 


SUMMARY 


uncoll. 

TOTAL 

HUILOI'R 

PaCKMARD 

FORHARD 

T|$*Ut 

nose  (RADS*CM*»? 

) 

3,91*09 

5.08*12 

1.30*03 

1.85*12 

3.92*12 

SlllCOK 

DOSE  <RA0S*CM»*2 

) 

9,01*08 

9,80*10 

2.39*0? 

1 .99* 1 0 

7.82*10 

FLUtNCFl  N 

) 

9.05*17 

8,22*21 

9.08*03 

3.50*21 

9.81*21 

CURRENT!  H 

) 

9,05*17 

7.55*20 

8.35*02 

-1 .72*21 

2.98*21 

f\t  ur.y 

FLU*  NCE (  N-MpV 

) 

8.19*18 

3.15*21 

5.13*0? 

8.20*20 

2.33*21 

hmi.v 

CURRENT (  N-MfV 

) 

8,19*18 

1,17*21 

1.91*0? 

-3.79*20 

1.55*21 

EXT.  E  L  X ,  «T.  (RADS 

1 

9.05*17 

8.15*21 

9.01*03 

3.59*21 

9,81*21 

AVERAGE 

ENERGY(M£V  1 

8,79*00 

3,89-01 

5.85-02 

2.31-01 

5.88-01 

NllitHtlN  GAMMA  SUMMARY 


iinlolL. 

total 

Ill'll  III  ip 

RACknaRO 

EORMARI' 

ITSSI'l  POSE  (9A0S*CM**2 

)  0.00 

5.50*12 

0,00 

8.82*  U 

5.05.12 

S1LKDN  OPSE  <HAnS«CM**2 

)  0.00 

8,03*12 

0,  JO 

7.91*11 

5.29*12 

MJMUER  ELUENCFl  G 

)  o.on 

1,39*22 

0.00 

5.5  3*21 

9.98*21 

NUMHl'R  CuH«rNT(  r. 

)  0.00 

3,83*21 

0.00 

•2.66*21 

6.29*21 

E^Rl'.Y  FLIlkNCEt  G-*EV 

)  0,00 

1 .55*2? 

0,00 

1.09*21 

1.99*2? 

ENIMGY  CURRENT!  G-mEV 

>  o.oo 

1 .23*22 

0,00 

•9,59*20 

1.27V22 

EXT,  ll  X,  *1,  ( 

)  o.oo 

o.oo 

0,00 

0,00 

o.oo 

AVERAGE  E NERGY ( I’EV  ) 

0,00 

1.11*00 

0.00 

1.87-01 

1 .95*00 

♦  ♦EXECUTION  COUPLE  TEO 

♦  STnP 
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SAMPLE  PROBLEM  2 
Page  1 


•  TlTtf  SAHPtf  AfW  PR0HLE**  2 

•  *j- souwr t  (  j) 

•  an  «s 
*ns.«“.  iu 

•  e»c 


ATM  PRORLEm  NUMRPR  I  SAMPLE  AT*  PRORLEM  ? 


NEUTRON  SOURCE  !  NTE  WNAL  l«**tV 

MOM "  AL 1 7AT I  Ol.c  1.000*00  NEUTRON  /KT*  YIELD*  l, 000*00  *T 
TOTAL  OUTPUTS  1,000*00  NEUTRON 

SOURCE  SPECTRUM 


ENERGY  (Rf  V  ) 

N/RT 

N/KT  *Mf  v 

1 .07-05- 

2.90-05 

0 ,  no 

0.00 

2. 90-05- 

1 .01-0* 

■  0.00 

0.00 

1  .01-09- 

5.85-09 

0.00 

0.00 

5.EE1-0*- 

1.15-01 

0.00 

0.00 

1.35-03- 

1.11-01 

n.oo 

0.00 

1.11-01- 

s.so-ot 

0.00 

0.00 

s.so-ni- 

1.11*00 

0.00 

0.00 

1,11*00. 

1 ,R3*no 

0.00 

0.00 

I.B1*C0- 

2.15*00 

0,00 

0.00 

ENFRGY(MEV  ) 

N/KT 

N/KT»ME  v 

2. 15*00. 

2,96*00 

n.oo 

o.co 

2.96*00. 

1.01*00 

n.oo 

n.oo 

1,01*00. 

a. o7*oo 

0.00 

0.00 

a, 07*oo- 

a, 97*00 

n.oo 

0.00 

9.97*00- 

6.16.00 

0.00 

0.00 

6,  16*00- 

B. 19.00 

0.00 

0.00 

6.19*00- 

1.00*01 

n.or 

0,00 

1.00*01- 

l .22*01 

o.oo 

0.00 

1 • 22*  0 1- 

1.50*01 

1.00*00 

1,57-01 

ATM  PROBLEM  MUMPER  1  SAMPLE  AT*  PROBLEM  2 


GHOUNO  LEVEL  . 00C* m ,  . 0 0OGm/C“* *2.  .OOOkET.  ,000«TLES 

*  mom  I E ,  RANGE.  RMS  7,071km.  62 . 1?9Gm/Cm* *2 •  21.199KET,  0.399“ILFS 

Slant  RANGE  MS*  lO.OOOKM,  52,989Gm/Cm**2.  32.606e>ET,  6,2iaMlLES 

•  TARGET  ALT,  NTs  27,071K“.1016,63uG‘VCM**2.  A0.6I6KET.  I6,B2tM?LES 

SOURCE.  ALT.  NS»  20.000RM,  979 . 1 69G-/C “« • 2*  65.617KFT.  12,«27Mll.E9 

SLAE.T  ANGLE  ANs  MN.OOOOEGRE  ES  (COS*  ,  7071  1  ) 

•  CALCUt  ATtO  PROM  OTMfR  COORDINATES 


ATR  PRo«Lf“  N1IMMER  1  Sample  AT*  PROBLEM  2 


NEUTRON  number  CURRENT 

DISTRIBUTION  OVER  ANGLE 

ANG  OtGREES  COSINE  S'JMQVER  ENERGY  CUEUHATIVF  SUM 


(  N/CM»*2 

•STER  )  t 

N2Cm»»2 

1 

160.0 

166,6 

•1,0000  - 

-.9728 

-7.129-15 

-1.299-15 

2 

166,6 

155.6 

-.9720  - 

-,9|00 

-7.072-15 

-9,016-15 

1 

155.6 

199.6 

-.9106  - 

-.815* 

•6,61 9- !  *■ 

-7,970-15 

a 

199.6 

1  11.7 

-.6159  - 

- . 6900 

-5.990-15 

-1.262-10 

5 

1  11.7 

122. B 

-.6990  - 

- .5*12 

-5,090-15 

-t .716-10 

6 

122.0 

111.0 

• , 59 1 ?  . 

-.1720 

-1.909-15 

•2.151-19 

7 

11  1 .8 

100. 9 

-.3720  - 

-.1899 

-2.521-15 

-2,990-10 

A 

100.9 

90,0 

-.1099  - 

.  DO  00 

•0.999-16 

-2.597-10 

R 

90.0 

n.  t 

,0600  - 

.  1  89* 

9,511-16 

-2,939-10 

10 

70,1 

60.2 

•  1 *99  . 

.1/20 

2,999.15 

-2.090-19 

1 1 

nR.2 

57.? 

,3720  - 

.5912 

5. 109-15 

-1.51 7-19 

1? 

57.2 

96.3 

.5912  - 

,69f'« 

8.919.15 

-7,258-15 

1  s 

96.1 

15.9 

.615* 

1. 195-19 

2.101-15 

1* 

15.9 

29.9 

,615*  . 

.9106 

1.707-19 

1.211-19 

15 

29,9 

11.9 

.9106  - 

.9720 

2.561-1* 

2.211-19 

16 

11.9 

.0 

,972*  - 

I ,0000 

6,796-19 

3,193-19 

90 


SAMPLE  PROBLEM  2 
Page  2 


{  N/C"»*i 


l 

0.00 

4.14-07 

7 

4, 1 4-0 7- 

1 .17-06 

1 

1 .17-06- 

}. flk-flk 

0 

5 . 0k-0k- 

1  . 07-oS 

s 

1.A7-AS- 

7. 01-OS 

k 

7.90-AS- 

1.01-04 

7 

1.01-04- 

S.M-04 

0 

S.k  1-04- 

J.lS-oJ 

9 

i.lS-01- 

1. ll-0| 

10 

1.11-01- 

S.SO-fll 

1 1 

s.sn-oi- 

1.1  MOO 

17 

1  • !  1  ♦ 0  3- 

1  .ki*oo 

11 

|  ,k  1*00- 

7, !S«00 

1« 

7.1S*00- 

7.  ak* 00 

IS 

7,46*00- 

1,01*00 

Ik 

i , o i ♦ oo- 

4.0  /«oo 

1  7 

4,0  7*00- 

4.97*0  C 

Ik 

4,97*00- 

k , Jk*  0  0 

1* 

k, |6*00- 

A. 19*10 

70 

0,19*00- 

1  .0,1*0  1 

71 

l .oo*oi- 

1.77*01 

77 

1.77*01- 

l.SO.fll 

Sum 


9.S74-I I 
«.74i-i i 
i  .Sh7-i i 
5 , 7*k-  )  f 

1 /0-i«» 
S.4P7-I < 
I  .0  1V-1  1 
1  .  '•  7  1 4 

4  ,Ok9-lS 
4.SM-1S 
?  .  H  1  «-  1  •, 
«,Vb4-1S 
7.7S7-1S 
4 .710- IS 

I . i<*a-is 
7.  St  l«|S 
I  .4S7-1S 
1.1  I4-IS 
S.939-|k 
k,*»»l-IS 

i.na-is 


0  I  S 1  w  1  H  i>  T  II)-,  iwtk  rululiv 
)  S  1J“  OVIM  4M.il  CIIMUL » I  I  VI 

•“In  )  (  II  /C  M  •  ♦  7 

,  1m  1- 1 0 


S. /7M-I  / 
I  .«»«*»- 1 fc 
2. 060-16 
1.761-16 
4,170-16 
S./lH-|k 
0.1H7-I6 
I.I77-IS 
7.777-IS 
4,077-IS 

k. k?7-is 
D.6SI-IS 

l.  0/1-14 

I.1S7-M 
I..1M9-I4 
1  .SI  7-1  4 
I .  /u-n-ia 
7.01 7-14 
7.774-14 

7.S70-I4 


••mcull  0*1  CPMMLfTH) 

•  link 
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SAMPLE  PROBLEM  3 
Page  1 


•  nut  s»m»'u  ath  PHOHitM  i 

•  •■-SOOmCF  (9) 

•  i-sv»i  noon  .015s  .tm  ,198  ,2m  .1?  .0/1  .009  .107  ,o!5  ,015  ,oi« 

,OM  .00089  ,00016 

«>-s  iooo 
♦Wh  1S00 

•  "TtP'».  1  9  5  7  to  IS  20 

•  oust  /m;/ 

•  f  *Ct 


ATR  PROBLEM  NU“SER  t  SIMPLE  ATP  PROHLtM  J 


HCUTRCn  SOURCE  EXTERNAL  HISTOGRAM 
NPRHALlZATjOua  1,002*00  NEUTRON  /KT  *  TIELDa  1,090*00  KT 
TOTAL  OUTPUT*  1.002*00  NEUTRON  /KT.  Y! 


SOURCE  SPFCTftUM 


ENERGY tMEv  ) 

N/K  t 

N/KT««EV 

f NE RGV ( me v  ) 

N/«T  1 

V/KT*Mf  V 

1 .07-05- 

/.90-OS 

0.00 

0.00 

2.15*00- 

2,96*00 

2.10.0/ 

2.09-01 

/ .90- C  S- 

1.01-09 

n.oo 

0.00 

2,96*09- 

1.01*00 

8.90-02 

1.62-01 

1,01-00- 

5.61-09 

0.00 

0.00 

1,01*00- 

9.07*00 

1.07-01 

1 .01-01 

5.81-fla- 

1.1S-01 

0.00 

0.00 

9.07*00- 

9 , 9  7  ♦  0  0 

1.50-0/ 

1.89.0/ 

1.15-01- 

1.11-01 

1 . 55-02 

1.99-01 

9.97*00- 

6.16*00 

1.50-0/ 

/.5/-02 

1.11-01- 

S.SO-01 

1.91-01 

1.26-01 

6,16*00- 

8. 1 9* On 

1.90-0? 

7,65-01 

5.5P-P1- 

1.11*00 

1,90-01 

1. 59-01 

«.|9*00- 

1 .00*01 

1.00-01 

1 .66-01 

1 . 11*00- 

1 ,61*00 

2.10-01 

I.Ol-oi 

1.00*01* 

1  ./2*01 

A.Oo.Ou 

9,05-OU 

1  ,M.I)I>. 

/.  15*00 

1 .20-01 

2.11-01 

1  ,??*01- 

1 .50*01 

1 .60-09 

5.71-05 

A  TR  P»OPLFm  num»ER  1  SAMPLE  ATR  PROBLEM  1 


ground  l*vei 

mOOI/,  RANGl  OM« 

•  Slant  range  rs* 
t  a  hue T  alt,  MTa 

SOURCE  AIT.  HS» 

•  slant  angle  am 


,000«m,  ,OOOGm/Cm»*2.  .OOORFT, 
1.500KM,  v’o, J71G-/CM»*2,  R.921KFT, 

l.iOOKM,  llh.  17lG-/CM»»2.  9.9/1KFT, 

1,0C"«M,  }|8.<j6/Gm/Ch**2.  9,8o/*U. 

1,000km,  1I8.962Gm/Cm»*2.  9.892KFT, 

.COODEGRECS  (CUSa  1,00000) 


•CALCULATED  FROM  OThIR  coordinates 


.OOOMJLES 
. 9  3/m  j  lF  5 
, 9  1/m i lF  S 
1 .»6«“!LF  S 
1 ,  OMim  ILE  S 


NEUTRON  Gamma  ROSE  VS.  target  alt,  (RADS*Cm»»?  ) 

HOSE  is  HE  NOE  R  SON  TISSUE 

OPSE  /■  CONCRETE 

OOSE  la  AIR 

OOSF  91  IONIZING  SILICON 


TARGET  ALT, 
1,000*00  KM 
9.000*00  km 
S. 000*00  KM 
7,000,00  km 
1,000.01  KM 
I  . S  o  0 ♦ o  t  KM 
/, 000*01  »  M 


OOSF  J 
i ,M 7-  l  1 
1 ,/T-ll 

H. 0/-1/ 
S.  i  ll-  1  1 
i./n-m 

S.Oh-  1* 

I. 17-1* 


POST  2 
/.OS-11 
1.3S-I1 

<i  .«*- 1  / 

T.3/-1  1 
I , 1/- 1 9 
N. 58-16 
1 ,/H- 1 0 


DOST  1 
l ,99-1  I 

I .os-i  1 

1.5/-12 
/.0/-1  1 
1.00-19 
9  .  ■./.  1  6 
1 ."2-16 


OOSE  9 
2,00-t  1 
1.  12-11 
9.11-1/ 
1.22-1  l 
1 ,?S-19 
5.91-10 
I .25-|h 


•  *t  *1  CUT  10n  CO-PCF  TCP 
•MT  1000 

•Rm.km,  ,s  |  |,s  22, S  1  «  S 

♦  t  »C. 


1  I 


I 

SAMPLE  PROBLEM  3 
Page  2 


ATH  PMORLfcX  NL'MPfR  2 


SAMPLE  ATR  RROMltM  3 


ground  iivTl 

Mfiuw.  NX* 

•  si  a*.t  #ahcf  w a» 

TAUf.FT  AIT,  HT« 
3  O'  i*-  C  A  ALT.  i"i« 

•  Sl«hT  ANGLF  A>)« 


.OflOKH.  ,OOOGm/c>«**2*  , OOOKF  T » 

.SOOt**  AS.aS70n/Cm*»2,  t.PAnKM, 
.SOOKM,  an. us  (G“/C“**2»  1.6#0KMi 

S.flPOK".  318.A62GM/cm*»2*  M^KMi 
J.nCOAx,  3tA.M62GM/CM»*2»  A.BA/xM* 

,(,oon> gaffs  (cos*  1  .ooooo) 


•  c*icui  *ifo  from  othfr  codrdinatfs 


.P06MTI  ES 

.3IIMRFS 
.311X11  fs 
t  .PMmll  f  S 

1.86am  1 1  FS 


NfllTRGN  GAM«A  DOSE  VS,  HOM12.  MANGE  (RA0S*CH»»i  ) 

oust  i»  xtNiHHS'm  tissue 

onst.  2»  Ct'*»CWP  T  r 

DOST  1«  AIM 

ousr  «■  IONIZING  SILICON 


HOXIZ.  RANGE 

OOSf  I 

DOSt  2 

OOSF  3 

OOSF  A 

5.000-01 

HN 

5,60-1  1 

5.RS-M 

#.96-11 

5,89-11 

1.000*00 

** 

A. Si-1 1 

A.R2-11 

A. 01-11 

A. 83-11 

1.500*00 

K*1 

1 .87*11 

2.05-lt 

1.60-11 

2.00-11 

2.000*u0 

A* 

6.  PA-12 

7.61-12 

6.02-12 

7.39-12 

2.500*00 

KM 

2.A2-I2 

2.69-12 

2.12-12 

2.61-12 

3.000*00 

¥  M 

8.AO-13 

9.32-13 

7.3A-13 

9.03-13 

A, 000*00 

KH 

9.99-1  A 

1. 1 1-13 

8,73-1  A 

1.07-13 

5,000*00 

K  ^ 

I.2A-1# 

1 . 37- 1  a 

1.08-1# 

1 .53-1 A 

••EXECUTION  COXRLftEO 
•  STOP 


I 
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SAMPLE  PROBLEM  4 
Page  1 


♦  TITLE  SAMPLE  AT9  problem  4 

•M-SCURCl (i) 

*k'.s09h  l,n*22 

•  kj.Y t E t  0  100 

♦  C’S*KH.10 

•  CAN  -90(t0)90 
*U!NSTHA1NT/I./(RS  5  250  0) 

•£*C  ! 


ATM  PR(.«l.fc«  SUMMER  1  Sample  atm  PH()HLEiA  m 


(.IUTWOn  SOURCE  1NTFHNAL  1  MMf  v 

NOMMAL17AT  X  0  A- *  1.000*22  (F.PlMrr*  /AT*  YtELO*  1.000*02  KT 
TOTAL  (IUTPUT*  l,flfl0*2ii  nEuTw(i*i 

source  spectrum 


FNF  MGF  {  *'E  V  ) 

N/«  1 

i  /At*  *1  V 

F  (<FMC.Y  (('EY  ) 

N/KT 

N/KT*MfV 

1 ,«7-(lS- 

2,90-05 

0.00 

0,0  0 

2.15*00- 

2.40*00 

0,00 

0.00 

2.9o«ns- 

1. 01-00 

o.oo 

0.00 

2 , MO*  0  0- 

1.01*00 

0.00 

0,60 

1,01-04- 

5,0  3-04 

0.00 

o.on 

1.01*00- 

4,07*00 

0.00 

0.00 

S.Ml-OU- 

■  *.15-01 

0.00 

o.oo 

4.07*00- 

•1.47*011 

0.00 

0.00 

3.15-ni- 

1.11-01 

O.on 

ii  ,oo 

4.97*00- 

A  .  30*  (l« 

0.00 

0.00 

1.1 1-01- 

5.50-01 

0.00 

o.oo 

o , Is*  oo- 

0. 19*i)0 

0,00 

0,00 

5.S0-0  1  - 

1.11*00 

o.oo 

0.00 

M. 14*00- 

1.00*01 

0,00 

0,00 

1.1 1*10- 

1.03*00 

0.04 

0.00  • 

1.00*01- 

1.22*01 

0,00 

0.00 

1  .  *1  1*00- 

2.15*00 

O.on 

0.00 

1.22*01- 

1.50*01 

1.00*24 

1.57*23 

,,,,,,  ,,,,,,, 

ATM  PRUMLE*-  MJMMtM  1  SAMPLE  atm  PNflPLM  4 


M  i'TROi 

60 STRAIN! 

'  FUR 

SLANT  RANGE 

A  NO  HFSi'IMSON  TISSUE*  2,50*01 

♦mO!7,  MA'.f.t 

•  Slant  mai6,f 

♦  T  A  rtf*  M  ALT, 

SOURCE  alt. 

SLANT  AnGLF 

(.,15  7  n-05 

1  ,4904*0  3 

¥ 

«.nons*ii  3 

1.0000*04 

-90.000 

DFG 

1,<»7  i  n  ♦  n  2 

1 .99(3  3*0  3 

V 

H.O  32t*C.  3 

1,  '000*04 

-SO, 000 

dec; 

a,m2a'**02 

1 ,9959*0  3 

► 

c.  1245*03 

1  ,  '1000*04 

-70,000 

OF  G 

1.0670*03 

2.01  M(,«  111 

M 

H.255S*  ii  3 

1  ,  0000*04 

-60.000 

DFG 

1  .  *  0  3  n  ♦  ij  3 

2,62n*0  3 

M 

H,  44  7  1  *  0  3 

1.0000*04 

-50.000 

DE  G 

1.56-0*03 

2.0 4 'ip* •!  3 

1 

H,hilA0*03 

1  ..1000*04 

-Mil,  ooo 

Off, 

1  , /M<2*  03 

2. "A 4r*n 3 

M 

a,  . n  7a* o  3 

1.0000*04 

-30,000 

DEG 

1 .4*.  1  1  *0  3 

|2 , OH A9 ♦  r  3 

¥ 

4,2HA2*nl 

1.0000*04 

-20,000 

OF  G 

2.07HM.O  1 

2, 1  1  09*1’  3 

¥ 

4. A  31u*0 1 

1,0000*04 

•  10,1104 

Of  G 

2,1400*03 

2, 1  40  0*11 1 

¥ 

1  .0110  0*0  4 

1 . 0000*04 

,000 

DEG 

2.11  36*6  1 

2, 1  A  A  5  ♦  3 

► 

1,03/6*04 

1,0000*04 

lo.ooo 

DEG 

2, OP n 5*03 

2.1427*03 

¥ 

1 .0750*04 

l ,  4000*04 

20.004 

ore 

1  .  5  2  ii  *■  ♦  0  3 

?.?1  7*.*6  3 

M 

1  , 1 109*04 

1  ,oono*04 

10,000 

OEC 

1  .711.2*01 

2,2404*01 

M 

1,1440*04 

1,0000*04 

40.000 

DFG 

1.4529*01 

'2,2602*0  3 

¥ 

1,1731*04 

1.0000*04 

50,000 

PFG 

1,1171*03 

2.2/4A*m 

¥ 

1,1470*04 

1  ,0000*04 

60,000 

OEG 

7,0190*02 

2,2420*03 

H 

1,2154*04 

1,0000*04 

70,000 

DEC 

4.0011*02 

2,3042*03 

¥ 

1  ,2269*04 

1,0000*04 

SO, 000 

OEG 

7.1  019-05  2,  3009*01  M 

♦  CALCUI  ATEO  EHdlM  OTHER 

1,2  307*04 
COORDINATES 

1 ,0000*04 

90,000 

DEG 

♦♦EXECUTION  COMPLETED 
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APPENDIX  A 

A  bibliography  of  basic  data  for  radiation  transport  in  air  has 
been  compiled  during  the  course  of  this  work  from  published  and  unpub¬ 
lished  works  which  have  resulted  during  the  past  fifteen  years.  This 
compilation  represents  the  major  comprehensive  works  of  which  the 
authors  are  aware,  based  on  a  systematic  document  search  by  computer, 
personal  communications  and  private  files. 

A  summary  of  references  on  theoretical  and  experimental  work 
is  listed  in  Tables  I  and  II,  respectively,  in  chronological  order,  accord¬ 
ing  to  type  of  radiation  considered.  Tables  III-V  give  details  of  the  calcu- 
lational  problems  considered  for  exponential  air,  infinite  uniform  air  and 
for  an  air-ground  interface,  respectively. 
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Table  L  Summary  of  References  on  Calculations!  Air-Transport  Data 
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1970  Striker'0'  Straker 


Table  IV.  list  of  Calculations  for  Infinite  >\lr 
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*>  denotes  fluoncc,  g  denotes  secondary  gamma  ray,  D  denotes  dose. 
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Table  IV.  (continued) 
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50 

d(r,E) 

D(r) 

SORS 

V  denotes  fluence,  g  denotes  secondary  gamma  ray,  D  denotes  dose. 
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Table  V.  List  of  Calculations  for  Air/Ground 


Investigator 

Date 

Method 

Source 

Energy 

Source 
He'L  ht 

Pages 

Results* 

Code 

Bigcers, 

Brovrn,  Kohr4U 

i960 

Monte 

Carlo 

12  Energies 
Weapons  Leakage 

300  ft 

173 

«(r,E,t)+ 

NHM 

Marcum42 

I960 

Of 

3,  lU  MeV 

300  ft 

55 

$(r*E) 

RAND 

Kinney 48 

1962 

Monte 

Carlo 

6  Energies 

0 

17 

d(r,E) 

05R 

Ritchie, 
Anderson4'*  »6b 

1962 

Monte 

Carlo 

1*»  Energies 

Iff— 

18 

♦(r,E) 

D(r,n) 

NHM 

Marcum  4 2 

1963 

Monte 

Carlo 

It*  MeV 

0,  750. 
3,000 

1i0 

«(r,E) 

D(r) 

CAPS 

DeVries61' 

1967 

Monte 

Carlo 

111  MeV 

116  ft 

75 

dfr.E.n) 

dg(r,E,fl) 

COHORT 

Straker  and 
Hynatf7* 

1967 

Discrete 

Ordinates 

Fission 

300  ft 

19 

❖( r ,E) 

D{r) 

DOT 

Straker 72 

1968 

Monte 

Carlo 

9  Energy 

Bands 

50  ft 

250 

d(r,E,n,t) 

D(r,t) 

0g(  r  ,E  ,0  ,t ) 
Dg(r,t) 

06R 

Straker72 

1968 

Discrete 

Ordinates 

9  Energy 

Bands 

50  ft 

250 

d(r,E,0) 

D(r,n) 

dg(r,E,n) 

DK(r,E,n) 

DOT 

d  denotes  fluence,  ft  denotes  secondary  gamna  ray,  D  denotes  dose. 


*Shock  vave  included. 
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APPENDIX  B 


DENSITY-SCALING  SOLUTIONS  TO  THE 
BOLTZMANN  TRANSPORT  EQUATION. 


Given  an  arbitrary  scattering  medium  and  particle  source  con¬ 
figuration,  the  corresponding  transport  problem  solution  can  be  obtained 
from  an  equivalent  (in  that  the  solutions  to  Boltzmann’s  equation  are 
identical)  density-scaled  configuration.  Relationships  are  given  here 
between  the  solution  to  the  given  arbitrary  configuration  and  the  solution  to 
the  transformed  configuration. 

The  Boltzmann  transport  equation  can  be  written 

3N(E,  «,?,  t)  ... 

- -t - =  S(E,  O,  r,  t)  -  v(E,  o)  .  7N(E,  f\  r,  t) 

-  at(E,  r)  •  p(r)  •  N(E,  P,  r,  t)  •  1  v(E)  ! 

4  ('  ]■  N(E\  C? r,  t)  ap(E ',"?)•  p(r) .  |v(E ')  I 
O’E' 

■  G(E\  P',  E,  °)  dE  dG  ' 

This  equates  the  time  rate  of  change  of  the  particle  differential 
number  density,  N(E,  0,r,  t),  to  the  source  strength  minus  the  losses  due 
to  convection  and  scattering  plus  the  gain  due  to  scattering  from  other 
energy  and  dii  ection  groups. 


o 

N  is-  the  number  of  particles  per  meter  -Mev-steradj.an 
S  is  the  source  of  particles  per  meter  -Mev-steradian 
is  the  total  particle  cross  section  in  meter  /target 

2 

o  is  the  particle-producing  cross  section  in  meter  /target 
P .  .  3 

p  is  the  target  density  in  targets /meter 

v  is  the  particle  velocity  in  meters/second 

G  is  the  probability  that  given  a  collision,  the  particle 
scatters  from  E',  ft'  to  E,  ft . 

If  we  perform  the  transformation 

Pt(?t)  =  p(r)/K 

rt  =  K  r 

tt  =  K  t 


Nt(E, 

6, 

?t,  v 

.  K-3 

N(E, 

$5, 

r,  t) 

St(E, 

ft, 

V  tt) 

=  K-" 

S  (  E , 

a. 

r,  t) 

then  the  transport  equation  is  merely  multiplied  on  both  sides 
-4 

by  K  .  In  other  words,  if  N(E,  ft,  r,  t)  is  a  solution  of  Boltz¬ 
mann  equation  given  p(r),  S(E,  ft,  r,  t),  then  Nt(E,  ft,  ,  tt) 
is  a  solution  given  the  configuration  pt  and  St(E,  ft,  ?t,  t.) 
as  defined  in  the  above  transformation. 

Information  about  the  particle  density  o  ( r,  t)  in  the  initial  problem 
is  equivalent  to  information  about  the  particle  density  it  a  scaled  time  and 
position  in  the  transformed  (scaled)  problem.  The  scaling  of  specific 
quantities  derived  from  the  number  density  follows  from  their  definitions. 
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